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demand,  and  catalyst  degradation.  Additionally,  an  empirical  model  is  developed  that 
predicts  the  fuel  conversion  in  the  reactor  for  a  given  set  of  fuel  flow  rate,  catalyst  bed 
length,  and  acoustic  pressure  parameters. 

Although  the  test  facility  used  was  not  optimized  for  utilizing  acoustic  waves, 
significant  acoustic  enhancement  of  the  steam-reformation  process  is  demonstrated  by  the 
study  and  such  enhancement  is  statistically  proven  to  have  a  positive  effect  on  the 
reformation  process.    Potential  benefits  resulting  from  acoustic  enhancement  of  steam 
reforming  as  shown  by  this  study  are  an  increased  reactor  capacity  for  a  given  size  and 
mass,  smoothing  of  the  temperature  profile,  better  temperature  control  of  the  catalyst  bed, 
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Conclusions  and  recommendations  for  further  study  and  application  are  drawn  from  the 
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CHAPTER  1 
INTRODUCTION 

This  chapter  gives  the  background  for  the  present  study  and  is  an  introduction  to 
fuel  cells,  including  a  brief  historical  sketch  of  the  technology.  The  different  types  of  fuel 
cells  are  outlined  and  the  main  difficulties  in  applying  fuel  cells  for  transportation  are 
mentioned.  Fueling  options  for  fuel  cell  vehicles  are  discussed  and  the  difficulties  in 
reforming  liquid  fuels  are  presented.  This  background  allows  a  formulation  of  the 
problem  definition  and  a  statement  of  the  study's  objective. 

Background 

Because  of  high  energy  conversion  efficiencies  and  low  emissions,  fuel  cell  power 
plants  (or  "engines")  are  currently  receiving  a  high  degree  of  interest  as  potential 
replacements  for  the  internal  combustion  engines  in  transportation  applications. 
Fuel  Cells 

The  fuel  cell  was  first  demonstrated  in  1839  by  Sir  William  Grove.    However,  the 
majority  of  practical  work  on  fuel  cells  did  not  happen  until  Francis  Bacon  demonstrated 
a  5  kW  working  system  in  1959.  NASA  became  interested  in  the  technology,  and  fuel 
cells  were  further  developed  under  General  Electric  and  United  Technologies  Corporation 
for  the  Gemini  and  Apollo  space  flights.  The  Arab  oil  embargo  in  1973  stimulated  the 
United  States  Department  of  Energy  (DOE)  to  look  for  ways  to  conserve 
petroleum-derived  fuels  and  utilize  fuel  cells  in  terrestrial  applications.    A  major 
breakthrough  in  fuel  cell  technology  came  in  the  1980s  when  Los  Alamos  National 

1 


Laboratory  and  Texas  A&M  University  found  a  way  to  decrease  the  platinum  loading  for 
certain  types  of  fuel  cells  (Norbeck,  1996).  Since  that  time  vehicle  manufacturers  have 
been  looking  to  the  fuel  cell  power  plant  as  a  way  to  comply  with  ever-increasing 
emissions  standards,  yet  keep  the  high  level  of  performance  expected  with  the  mature 
internal  combustion  engine. 

The  fuel  cell  is  an  electrochemical  cell  that  uses  an  oxidation  reaction  to  generate 
electricity  directly  from  a  fuel  and  an  oxidizer.  Figure  1 . 1  shows  how  an  acid  or  polymer 
electrolyte  fuel  cell  operates.  Other  fuel  cell  types  operate  in  a  comparable  manner  yet 
these  other  types  of  fuel  cells  can  use  conductive  oxide,  alkali,  or  molten  carbonate  as 
electrolytes.  For  the  acid  or  polymer  type  fuel  cells  as  shown  in  Figure  1.1,  hydrogen 
fuel  is  brought  to  the  anode  of  the  fuel  cell.  At  the  anode,  typically  a  platinum-type 
catalyst  splits  the  hydrogen  into  protons  and  electrons.  The  protons  are  free  to  migrate 
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Figure  1.1:  Simplified  fuel  cell  schematic. 


across  the  electrolyte  and  combine  with  oxygen  at  the  cathode  to  form  water.  To 
complete  this  reaction,  the  electrons  that  are  left  at  the  anode  must  be  brought  to  the 
reaction  site  at  the  cathode.  A  conductive  path  is  made  available  through  a  load  and  the 
electrons  are  passed  through  the  circuit.  The  end  products  are  water,  heat,  and  direct 
current  electrical  power.  Because  a  single  fuel  cell  only  generates  about  0.75  volts, 
practical  applications  arrange  multiple  fuel  cells  in  series  to  form  a  fuel  cell  stack  that 
yields  usable  voltages. 

Several  different  types  of  fuel  cells  exist.  The  different  types  of  fuel  cells  are 
normally  distinguished  by  their  electrolyte.  Table  1.1  (adapted  from  Norbeck.  1996)  lists 
several  types  of  fuel  cells  along  with  their  operating  temperature,  typical  application,  and 
chemical  tolerance  to  carbon  monoxide  and  carbon  dioxide. 


Table  1.1:  Types 

of  fuel  cells  and  their  operating  parameters. 

Fuel  cell 
type 

Typical  stack 

operating 
temperature  C° 

Typical  applications 

CO  tolerance 

C02 

tolerant? 

Solid  oxide 
(SOFC) 

1000 

Steam  cogenerating  power 
plants 

Good 

Yes 

Molten 
carbonate 
(MCFC) 

600 

Steam  cogenerating  power 
plants 

Good 

Yes 

Phosphoric 

acid 

(PAFC) 

1 50-205 

Smaller  stationary  power 

supplies 

Heavy  transport  (buses, 

trucks,  ships,  trains) 

Fair 

Yes 

Solid 

polymer 

(PEM) 

25-120 

Stationary  and  mobile  power 

supplies 

Heavy  transport,  light 

vehicles 

Poor 

Yes 

Alkali 

65-220 

Space  Vehicles,  submarines, 

specialty  vehicles  carrying 

hydrogen  and  oxidizer 

Poor 

No 

Vehicle  Applications 

Transportation  applications  place  demands  on  the  power  source  that  currently  rule 
out  certain  types  of  fuel  cells.  Some  of  the  more  significant  demands  on  the  fuel  cell 
system  for  these  applications  are  the  ability  to  use  air  as  the  oxidizer,  and  quick  start-up 
and  shutdown  periods.  Because  of  the  development  of  the  Apollo  and  Space  Shuttle 
Programs,  alkali-type  fuel  cells  have  been  thoroughly  studied,  produced,  and  successfully 
operated  for  over  40  years.  These  fuel  cells  have  provided  astronauts  with  electrical 
power,  heat  and  drinking  water.  Excellent  as  they  may  be  for  such  space  vehicles 
carrying  hydrogen  and  oxidizer,  these  types  of  fuel  cells  are  of  little  use  in  terrestrial 
applications  because  the  Alkali-type  fuel  cells  are  not  tolerant  of  the  low  levels  of  carbon 
dioxide  in  the  air  or  in  reformed  liquid  fuels.  Solid  Oxide  and  Molten  Carbonate  type 
cells  are  excellent  in  terms  of  carbon  dioxide  and  carbon  monoxide  tolerance;  however, 
their  high  temperatures  and  associated  long  start-up  times  currently  prohibit  their  use  in 
automotive  applications.  The  lower-temperature  fuel  cells  that  are  carbon  dioxide 
tolerant  such  as  the  Polymer  Electrolyte  Membrane  (PEM)  and  Phosphoric  Acid  Fuel 
Cells  (PAFC)  are  more  suitable  for  transportation  applications.  PEM  (then  referred  to  as 
solid  polymer  electrolyte  or  SPE)  fuel  cells  were  first  developed  in  the  1 960s  for  NASA 
and  were  used  in  the  Gemini  Program  (Prater,  1990).  These  fuel  cells  were  not  further 
developed  for  NASA  as  alkali-type  cells  were  used  in  later  space  missions.  PEM  fuel 
cells  can  have  high  power  densities  because  of  the  thin  polymer  electrolyte  and  thus  have 
an  advantage  over  the  bulkier  PAFCs  for  terrestrial  applications.  On  the  other  hand, 
PAFCs  which  operate  at  a  higher  temperature  have  a  better  tolerance  for  carbon 
monoxide  in  the  fuel  stream.  Both  of  these  cell  types  (indeed  all  fuel  cells,  except  for 


direct  methanol)  best  utilize  pure  hydrogen  as  the  fuel.  As  it  has  little  infrastructure,  this 
use  of  hydrogen  is  a  significant  drawback  to  widespread  implementation  of  hydrogen- 
fueled  fuel  cell  vehicles.  As  shown  in  Figure  1.2  (Photograph  by  Tom  LaRocque  and 


Figure  1.2:  Hondas"  hydrogen-fueled  experimental  fuel  cell  vehicle. 
NREL/DOE),  this  limited  availability  has  not  prevented  development  and  demonstration 
of  hydrogen  prototype  fuel  cell  vehicles  by  most  auto  manufacturers.  Auto  manufacturers 
see  fuel  cell  powerplants  as  perhaps  the  only  feasible  technology  that  can  match  the 
requirements  of  both  consumers  and  increasingly  stringent  emission  regulations  while 
reducing  dependence  on  petroleum  and  reducing  the  generation  of  C02. 
Fueling  Options 

In  a  fuel  cell  vehicle  that  uses  hydrogen  as  a  primary  fuel,  there  are  no  vehicle 
emissions  as  the  only  end  products  are  water,  heat,  and  electrical  power.  As  hydrogen  is 
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also  the  most  abundant  element  in  the  known  universe,  it  may  seem  that  hydrogen  would 
be  the  ideal  primary  fuel  for  a  fuel  cell  vehicle.  Unfortunately,  virtually  all  of  the 
hydrogen  on  earth  is  bound  up  in  compounds  such  as  water  that  are  at  a  low  energy  state. 
Because  there  are  no  direct  hydrogen  sources  available,  with  the  exception  of  small 
quantities  of  hydrogen  found  with  some  natural  gas  deposits,  hydrogen  fuel  must  be 
derived  from  other  sources.  Although  there  are  some  projects  studying  and 
demonstrating  renewable  hydrogen  production  (Erickson  and  Goswami,  2001),  the 
majority  of  hydrogen  production  currently  utilizes  fossil  fuels,  primarily  natural  gas,  for 
the  feedstock.  The  production  of  hydrogen  from  fossil  fuels  creates  emissions  where  the 
hydrogen  or  the  energy  to  produce  the  hydrogen  is  generated.  For  a  hydrogen  fuel  cell 
vehicle  there  are  no  vehicle  emissions  but  the  production  of  the  hydrogen  fuel  creates 
emissions  upstream  for  the  vehicle.  While  there  are  places  such  as  the  Los  Angeles 
Basin,  where  transferring  pollution  outside  of  the  area  would  be  beneficial,  the  overall 
efficiency  is  changed  only  slightly  by  production  of  hydrogen  from  a  hydrocarbon  fuel 
on-board  the  vehicle  (Louis,  2001).  Adamson  and  Pearson  (2000)  have  shown  that  the 
"well-to-wheel"  efficiency  and  emissions  for  fossil-derived  hydrogen-fueled  fuel  cell 
vehicles  over  the  entire  use-phase  life-cycle  are  only  slightly  better  than  that  of  fossil- 
derived  methanol-fueled  fuel  cell  vehicles  over  the  same  life  cycle. 

With  these  similarities  in  the  overall  energy  conversion  understood,  one  must 
consider  which  fuel  would  be  a  wise  choice  for  use  in  fuel  cell  vehicles.  Table  1 .2 
(Adapted  from  Ogden  and  Nitsch,  1993;  and  Turns,  1996)  shows  the  energy  densities, 
based  on  the  lower  heating  value  available  from  various  fuels.    These  energy  densities 
demonstrate  a  practical  area  where  hydrogen  fuel  is  clearly  not  the  best  choice.  While 


Table  1.2:  Energy  densities  of  fuels  (LHV). 

Fuel 

Volumetric  energy 
density  (M J/1  iter) 

Mass  energy  density 
(MJ/kg) 

Diesel 

33.3 

44.5 

Gasoline 

32.3 

43.0 

Liquified  propane  (690  kPa) 

23.6 

46.3 

Methane  gas 
(16.5  MPa) 

5.51 

50.0 

Liquified  methane 
(-161  °C) 

21.5 

39.1 

Methanol 

15.9 

19.9 

Ethanol 

21.6 

27.2 

Hydrogen  gas  (16.5  MPa) 

1.76 

120.2 

Liquified  hydrogen  (-253  °C) 

8.4 

120.2 

Metal  hydride 

4.9 

1.80 

hydrogen  gas  has  the  highest  specific  energy  (mass  energy  density)  of  all  practical  fuels, 
it  has  the  lowest  volumetric  energy  density.  The  volumetric  energy  density  of  hydrogen  is 
still  paltry  compared  to  liquid  hydrocarbon  fuels  even  in  cryogenic  liquid  form  near 
absolute  zero  temperature.  This  issue  of  poor  volumetric  energy  density  suggests 
prohibitively  large  fuel  tanks  for  hydrogen  fuel  cell  vehicles  to  achieve  current  vehicle 
ranges  (-560  km,  -350  miles).    In  terms  of  volumetric  energy  density,  a  more  feasible 
alternative  is  to  use  a  liquid  hydrocarbon  fuel  such  as  liquified  methane,  liquified 
propane,  methanol,  ethanol.  gasoline,  or  diesel. 

As  previously  mentioned,  an  even  more  significant  problem  hindering  the 
introduction  of  direct  hydrogen  fuel  cell  vehicles  is  the  lack  of  hydrogen  fueling 
infrastructure.  The  current  United  States  fuel  infrastructure  has  evolved  for  use  of  liquid 
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fuels,  notably  gasoline  and  diesel.  Because  of  the  infrastructure  and  energy  density 
issues,  near-term  true  consumer  (as  opposed  to  limited  fleet  operation)  fuel  cell  vehicles 
will  likely  utilize  a  liquid  hydrocarbon  fuel.  Although  smaller  environmental  benefits  are 
projected  for  using  these  liquid  fuels,  rather  than  using  hydrogen  from  a  renewable 
resource,  the  use  of  fossil-derived  liquid  hydrocarbon  fuels  can  serve  as  a  precursor  to 
enable  the  move  toward  a  renewable  energy  economy. 
Reforming 

Liquid  hydrocarbon  fuels  need  to  be  decomposed  or  "reformed"  into  a  hydrogen- 
rich  stream  before  they  can  be  used  by  a  Polymer  Electrolyte  Membrane  (PEM)'  or  a 
Phosphoric  Acid  Fuel  Cell  (PAFC).  This  can  be  accomplished  by  a  steam  reforming, 
partial  oxidation,  or  an  autothermal  reaction.  For  comparisons  of  the  reforming  methods, 
see  Ahmed  and  Krumpelt  (2001),  Brown  (2001),  and  Docter  and  Lamm  (1999).  The  on- 
board reforming  process  further  complicates  the  fuel  cell  vehicle  and  increases  the  cost, 
weight,  and  under-hood  volume  requirements  of  the  fuel  cell  engine.  The  reforming 
process  also  penalizes  the  overall  on-board  energy  conversion  efficiency  because  of  an 
associated  increase  in  entropy.    Furthermore,  the  liquid-fueled  vehicle  may  require 
clean-up  units  to  purify  the  reformate  stream  before  hydrogen-rich  mixture  can  be  utilized 
in  the  fuel  cell.  The  reforming  process  also  may  further  complicate  the  overall  vehicle 
system  by  requiring  long  start-up  (Erickson  et  al.,  2000)  and  shutdown  (Heckwolf  et  al., 
2001)  periods.  Although  liquid-fueled  fuel  cell  vehicles  are  demonstrated,  as  shown  in 
Figure  1 .3  (Photograph  by  Daniel  Betts  and  the  University  of  Florida  Fuel  Cell  Research 


1  There  are  some  exceptions  such  as  the  direct  methanol  PEM  fuel  cell  but  it  is 
still  in  the  laboratory  research  phase. 


Figure  1.3:  Methanol-fueled  PAFC  transit  bus  at  the  University  of  Florida, 
and  Training  Laboratory),  the  efficient  reforming  of  the  liquid  fuel  utilizing  practical  on- 
board reformers  is  considered  one  of  the  major  technological  hurdles  delaying  widespread 
introduction  of  fuel  cell  vehicles  (Erickson  and  Roan,  1999;  Romney.  2001 ).  Steam 
reforming,  where  fuel  and  steam  are  externally  heated  and  catalytically  reacted  is 
currently  the  most  widespread  method  of  reforming  a  hydrocarbon  fuel  into  a  hydrogen- 
rich  gas.  Although  steam  reforming  is  widely  used,  typical  steam  reformers  are  not 
designed  for  transportation  application  and  thus  they  have  large  volume  with 
corresponding  high  mass,  slow  transient  response  (especially  at  start-up),  and  are  often 
characterized  by  having  hot  and  cold  spots  throughout  (deWild  and  Verhaak.  2000). 

Problem  Definition 

As  introduced  above,  the  problem  of  reforming  liquid  hydrocarbon  fuels  to 
hydrogen  is  currently  one  of  the  major  technological  hurdles  delaying  the  introduction  of 
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liquid-fueled  fuel  cell  vehicles.  While  it  has  been  demonstrated  that  the  endothermic 
steam-reforming  process  can  be  accomplished  in  a  catalytic  converter  (Idem  and  Bakhshi, 
1994:  Ledjeff-Hey  et  al.,  1998:  Peppley  et  al.,  1997),  the  process  has  liabilities  of  weight, 
size,  and  complexity,  in  addition  to  long  start-up  and  transient  response  times  (Ohl  et  al., 
1996).  Because  of  these  liabilities  any  enhancement  scheme  for  the  steam-reforming 
process  which  could  lead  to  reducing  size,  weight,  or  response  time  of  the  reformer 
should  be  investigated. 

Research  Objectives 

A  steam-reformation  process  that  utilizes  acoustic  fields  in  critical  fluid  paths  is 
proposed.  This  study  develops  a  theoretical  approach  for  the  use  of  acoustic  fields  in 
conjunction  with  steam-reforming  of  methanol  and  demonstrates  that  the  resulting 
acoustic  enhancement  of  reforming  is  one  approach  that  can  reduce  the  above  liabilities 
of  the  steam-reforming  process. 

The  objectives  of  the  present  study  are  to  develop  theory-based  models  for 
acoustic  enhancement  of  the  reformation  process,  to  experimentally  demonstrate  a  steam- 
reformation  process  under  the  influence  of  an  acoustic  field,  and  to  quantify  the  effects  of 
acoustic  fields  on  the  steam-reforming  process.  This  includes  an  analysis  of  acoustic 
fields  and  investigation  of  the  effect  of  acoustic  fields  with  relation  to  the  reactor  output 
parameters  of  temperature  profile  and  controllability,  power  demand,  fuel  conversion, 
and  catalyst  degradation.  The  objectives  also  include  a  demonstration  of  the  overall 
hydrogen  power  output  relative  to  acoustic  power  input.  In  connection  with  the  objective 
to  quantify  the  effect  of  the  acoustic  field  on  the  reformation  process,  an  empirical  model 
is  developed  that  predicts  the  conversion  for  a  given  set  of  space  velocity,  and  acoustic 
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pressure  parameters.  Fulfillment  of  the  objectives  provides  understanding  of  how 
acoustic  fields  can  be  used  to  enhance  the  reformation  of  hydrocarbon  fuels,  resulting  in 
smaller,  more  effective  reformers,  especially  relevant  for  transportation  applications. 


CHAPTER  2 
LITERATURE  REVIEW  AND  THEORETICAL  APPROACH 

This  chapter  provides  a  theoretical  basis  and  background  for  the  strategy  of  using 
acoustic  waves  as  a  means  of  enhancing  performance  in  conjunction  with  the 
hydrocarbon  fuel  steam-reforming  process.  As  part  of  the  theoretical  approach,  the 
literature  dealing  with  steam-reforming  and  acoustic  enhancement  of  related  processes  is 
reviewed  and  the  potential  contribution  of  the  present  study  is  identified. 

Steam  Reforming 

Steam  reforming  is  the  conversion  of  a  hydrocarbon  and  steam  into  hydrogen-rich 
gas.  The  steam  reforming  process  was  developed  in  the  1920s  to  produce  a  mixture  of 
fuel  gas  containing  carbon  monoxide  and  hydrogen  known  as  "syngas"  or  synthetic  gas 
(Cox  and  Williamson,  1977).  Steam  reforming  can  occur  with  a  variety  of  hydrocarbon 
fuels  including  methanol  and  is  used  extensively  for  the  commercial  production  of 
hydrogen  from  natural  gas.  Equation  2.1  shows  the  ideal  steam-reforming  reaction  for 

CH3OH  +H20^>  3//2  +  C02  (2.1) 

production  of  hydrogen  and  carbon  dioxide  from  methanol  and  steam.    The  change  in 

enthalpy  of  the  reaction  is  +49.0  kJ/mol.  This  net  gain  in  enthalpy  indicates  that 

methanol-steam  reformation  is  endothermic.  Practical  product  species  typically  include. 

in  addition  to  hydrogen  and  carbon  dioxide,  carbon  monoxide,  unreacted  methanol  and 

steam.  Significant  methane  concentrations  and  solid  carbon  formation  have  also  been 

noted  under  certain  conditions.  The  reforming  process  typically  occurs  at  a  temperature  of 
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200-300  °C  and  in  the  presence  of  a  catalyst  (usually  copper  and  zinc  on  an  alumina 
substrate)  that  aids  the  reaction  kinetics  and  shirts  selectivity  to  carbon  dioxide.  Excess 
water  (1.3:1  -1.5:1  steamxarbon  ratio)  is  used  to  inhibit  formation  of  solid  carbon  (Cheng 
et  al.,  1999;  Amphlett  et  al.,  1993).    Models  have  been  developed  for  methanol  steam 
reforming  processors  (Asprey  et  al.,  1999)  and  reaction  mechanisms  have  been  proposed 
(Amphlett  et  al.,  1988).  It  has  been  proposed  that  a  simplified  model  of  the  methanol 
steam  reforming  process  can  be  seen  as  two  series  reactions  which  are  the  endothermic 
decomposition  of  methanol,  followed  by  the  exothermic  water-gas  shift  reaction.  These 
steps  are  shown  in  Equations  2.2  and  2.3  respectively.  The  decomposition  of  methanol  is 
thought  to  be  the  rate-limiting  step  as  its  rate  is  approximately  one  order  of  magnitude 
less  than  the  water-gas  shift  reaction  (Pour  et  al.,  1975). 

C//3O//0  CO+2H2  (2.2) 

CO+  H20<^>  C02  +  H2  (2.3) 

With  its  high  hydrogen  yield  from  a  liquid  fuel,  the  methanol-steam  reforming 
reaction  is  being  studied  as  a  possible  method  to  provide  hydrogen  for  fuel  cell  vehicles. 
Methanol-steam  reforming  has  been  shown  to  successfully  produce  hydrogen  for  vehicle 
applications  (Kaufman,  1992),  yet  the  process  is  hindered  in  transportation  applications 
by  characteristics  of  slow  transient  response  (Dusterwald,  1997).  temperature  gradients 
(Wild  and  Verhaak,  2000),  and  large  volume  and  mass  (Ohl  et  al.,1996). 

Limiting  Steps  in  the  Reformation  Process 
For  fuel  cell  vehicle  applications,  size,  weight,  and  transient  response  are  vital 
factors.  In  order  to  design  an  optimal  reformer  for  fuel  cell  vehicles  one  would  desire  the 
smallest  possible,  lightweight  reformer  that  would  have  enough  catalyst  surface  area  to 
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react  the  fuel  into  a  hydrogen-rich  stream  at  the  maximum  design  flow  at  the  end  of  the 
useful  lifetime.  Steam-reforming  of  hydrocarbons  can  be  limited  by  heat  transfer  rates, 
diffusion  or  mass  transport,  chemical  kinetics  and/or  by  degradation  of  the  catalyst  with 
time.  These  limitations  are  evidenced  by  restrictions  in  reformer  capacity,  slow  transient 
response  and/or  limited  lifetime. 

Because  of  its  endothermic  nature,  the  reformation  of  liquid  hydrocarbon  fuels 
requires  heat  transfer  to  the  reactants  and  the  catalyst.  In  a  steam  reformer  heat  comes 
from  an  external  source  and  is  transferred  through  the  catalyst  bed  housing  wall  and  into 
the  catalyst  bed  as  shown  in  Figure  2.1 .  Because  of  this  endothermic  nature  and  irregular 
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Figure  2.1:  Simplified  representation  of  a  packed 
bed  methanol-steam  reformer. 
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flow  patterns  associated  with  packed  bed  catalysts,  the  typical  steam  reformer  is  plagued 
with  temperature  gradients.  Typically  high  temperatures  are  found  nearest  to  the  heating 
surface,  and  low  temperatures  are  found  toward  the  midpoint  of  the  catalyst  bed. 
Furthermore,  if  the  heat  transfer  rate  to  the  midpoint  of  the  catalyst  bed  is  not  high 
enough,  unreacted  species  may  pass  through  the  reactor  at  this  location  resulting  in  a 
limited  reformer  capacity.  Additionally  the  catalyst  may  be  damaged  by  condensing 
species  where  heat  transfer  is  insufficient.  Heat  transfer  can  take  place  as  convection, 
conduction,  or  radiation.  Radiation  is  typically  negligible,  as  the  local  temperature 
difference  is  relatively  small.  Heat  conducts  through  the  catalyst  bed  housing  wall  and 
possibly  into  the  first  layer  of  catalyst  particles  coming  into  contact  with  the  wall. 
Because  of  point-to-point  contact  between  the  catalyst  pellets,  conduction  from  catalyst 
pellet  to  catalyst  pellet  is  minimal,  and  the  bulk  of  the  heat  transfer  to  the  center  of  the 
packed  bed  is  via  convection  through  the  fluid  inside  the  reactor.  As  the  explained 
temperature  gradient  manifests,  heat  transfer  via  convection  is  typically  the  limiting 
factor.  Heat  transfer  limitations  can  be  compensated  for  by  an  increased  size  of  the 
reformer  to  increase  the  heat  transfer  area  (Ohl  et  al.,  1996).  Such  design  results  in 
undesirable  large  processors  and  corresponding  high  mass.  The  large  size  and  thermal 
mass  then  penalize  the  process  for  transportation  applications. 

Diffusion  can  also  limit  the  reformation  process.  The  processes  of  transport 
within  the  catalyst  bed  and  conversion  to  a  hydrogen-rich  gas  is  represented  in  Figure  2.2. 
As  this  figure  shows,  the  hydrocarbon  and  steam  mixture  must  first  be  transported  to  the 
surface  of  the  catalyst.  Next  the  mixture  must  diffuse  through  the  pores  to  an  open 
catalyst  site  to  react.  The  reaction  rate  on  the  catalyst  is  then  governed  by  chemical 
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Step  1-  Reactants  diffuse  through  bulk 
stream  to  catalyst  surface 
Step  2-  Reactants  diffuse  through 
catalyst  pores  to  open  reaction  site 
Step  3-  Reactants  absorb  onto  catalyst 
Step  4-  Reaction  time  governed  by 
chemical  kinetics  on  the  catalyst 
Step  5-Products  adsorb  from  catalyst 
Step  6-Products  diffuse  through  catalyst 
pores  to  surface 

Step  7-  Products  diffuse  from  surface  to 
bulk  stream 


Steps  3-5 


Figure  2.2:  Simplified  representation  of  the  reactants  moving  to  the  open  catalyst  site 
and  products  moving  out  to  the  bulk  stream. 

kinetics.  Products  of  the  reaction  must  then  diffuse  back  out  of  the  catalyst  pores  and  into 
the  bulk  flow  stream.  Increasing  the  amount  of  catalyst  sites  available  can  compensate  for 
limitations  in  diffusion.  Normally  this  means  increasing  the  amount  of  catalyst  to  provide 
sufficient  reaction  area  and  open  reaction  sites  closer  to  the  surface  of  the  catalyst, 
allowing  for  a  decreased  diffusion  distance.  This  typically  requires  an  increase  in  the 
reactor  volume  or  conversely,  a  decrease  of  mass  flow  rate  of  the  reactants.  This  is  a 
decrease  in  what  is  known  as  "space  velocity"  as  shown  in  Equation  2.4.  This  decrease  of 
space  velocity  once  again  results  in  larger  reforming  processors  or  a  more  limited 
processing  capacity  for  a  set  reformer  volume. 
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Volumetric  Flow  Rate  ,_   .. 

Space  Velocity  =     ,,  , — (2.4) 

K  J  Volume  of  Reactor  v       ' 

Slow  chemical  kinetics  can  usually  be  alleviated  by  increasing  the  temperature. 
Because  temperature  gradients  exist  within  the  catalyst  bed,  temperature  cannot  be 
increased  at  the  midpoint  of  the  catalyst  bed  without  possibly  damaging  the  catalyst  near 
the  wall  surfaces  by  overheating.  Thus,  restriction  by  chemical  kinetics  is  normally  a 
manifestation  of  convective  heat  transfer  limitations  in  the  center  portion  of  the  catalyst 
bed. 

In  methanol-steam  reformation,  performance  decreases  as  the  catalyst  ages.  The 
lifetime  requirements  of  transportation  applications  require  the  reforming  processor  to  be 
over-designed  to  order  to  compensate  for  performance  degradation.  Once  again  this 
overdesign  results  in  larger  processors  and  increases  under-hood  volume  requirements  in 
transportation  applications.  Catalyst  lifetime  can  be  affected  by  poisoning,  fouling,  or 
sintering  of  the  catalyst.  Deactivation  due  to  poisoning  occurs  with  chemisorption  of 
impurities  on  reaction  sites  (Farrauto  and  Bartholomew,  1997).  In  reduced  copper-zinc 
catalysts,  oxygen  is  a  poison  as  the  metals  are  changed  to  an  oxidized  state  and  thus  are 
unavailable  for  the  steam  reforming  process.  Fouling  is  the  physical  blockage  of  reaction 
sites  and  can  take  place  when  fluid  phase  molecules  contact  the  catalyst  or  as  solid 
material  (e.g.,  carbon)  collects  on  the  surface  (Sterchi,  2001 ).  In  pelletized  catalysts 
where  the  interior  area  is  dominant,  fouling  at  the  catalyst  surface  or  pores  can 
significantly  change  the  available  catalyst  area.  Sintering  or  a  change  in  catalyst  structure 
normally  occurs  with  thermal  stress  of  the  substrate  material.  The  change  in  structure 
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reduces  the  effective  area  of  the  catalyst.  In  order  to  avoid  sintering,  overheating,  rapid 
thermal  cycling,  and  thermal  gradients  should  be  minimized. 

It  is  proposed  that  some  of  these  limiting  steps  in  the  methanol-steam  reformation 
process  may  be  minimized  by  imposing  an  acoustic  field  in  the  critical  flow  paths.  The 
proposed  benefits  due  to  this  acoustic  field  are  expected  as  a  result  of  increased  heat 
transfer  and  mass  transfer  rates. 

Acoustics  in  Past  Processes 

In  several  past  studies  (see  overview  by  White,  1991)  an  acoustic  wave 
superimposed  on  a  flow  field  has  been  shown  to  change  the  characteristics  of  the  flow 
field.  Periodic  plane  waves  in  a  wave  guide  can  be  represented  by  their  pressure  (p)  and 
velocity  (u)  components  as  shown  in  Equations  2.5  and  2.6  (Kinsler  et  al.,  1982).  In  these 
equations  x  is  position,  t  is  time,  j  is  the  imaginary  constant,  p0  is  the  non-oscillating 

(     „    2nf  \                 (          2*f  "| 
j  2*ft- x  j  2jift  + x  nSl 

A        j[27ift-  —  x)  B        j(2;ift+  — x) 

u  =  e  c    J  +  eV         c    }  (2.6) 

PoC  PoC 

density  and  A  and  B  are  found  from  the  initial  and  boundary  conditions.  A  general 
representation  of  an  acoustic  wave  at  a  frequency  (f)  in  combination  with  a  steady  flow 
component  within  a  tube  is  shown  in  Figure  2.3.  In  this  figure  the  acoustic  pressure 
(shown  as  the  shading  variation)  drives  the  acoustic  velocity.  When  the  acoustic  velocity 

component  (LT)  is  larger  than  a  left-to-right  steady  flow  component  (U),  the  motion  of  the 

fluid  periodically  changes  direction  depending  on  space  and  time.    The  use  of  resonance 
can  allow  for  high  acoustic  pressures  and  velocities  to  build  with  the  resulting  standing 
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Steady  Flow  (U)  with  a  Superimposed  Acoustic  Wave 
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Figure  2.3:  Flow  fields  at  different  times  for  a  plane  wave  superimposed  onto  a  steady 

flow  field. 

waves,  yielding  a  periodic  change  of  motion  although  now  dependant  on  the  spatial 

position  associated  with  the  specific  mode  of  the  standing  wave. 

In  the  development  of  fluid  mechanics  it  has  been  noted  that  acoustic  waves  can 

influence  several  parameters  of  the  flow  field.  These  parameters  include  the  boundary 

layer  and  flow  characterization.  Stokes'  2nd  problem  solved  for  the  velocity  profile  near  a 

flat  plate  with  an  oscillating  stream  (see  discussion  in  White,  1 99 1 ).  The  solution  of  this 

problem  allows  calculation  of  the  laminar  boundary  layer  thickness  in  this  oscillating 

flowfield.  The  boundary  layer  is  represented  by  Equation  2.7,  where  6  is  the  boundary 


8-6.5 


2*f 


(2.7) 


20 
layer  thickness,  v  is  the  kinematic  viscosity  of  the  fluid,  and  f  is  the  frequency  of  the 
wave  oscillation.  The  acoustic  field  may  also  change  the  flow  regime  as  discussed  by 
White  (1991 ).  The  kinetic  Reynolds  number  (Rek),  which  is  the  Reynolds  number  for 
oscillating  flowfields  is  shown  in  Hquation  2.8  for  pipe  flow.  In  this  equation,  r0  is  the 
radius  of  the  pipe  or  a  characteristic  length.  Like  the  standard  Reynolds  number  the  flow 
may  become  turbulent  when  this  value  approaches  2000.  Turbulent  flow  regimes  are 

Re,  =  ^^  (2.8) 

V 

characterized  by  chaotic  motion  of  the  fluid  resulting  in  increased  mixing,  and  for  pipe 
flow  a  flattened  velocity  profile  as  opposed  to  a  parabolic  profile  for  a  fully  developed 
laminar  regime. 

The  acoustically  thinned  boundary  layer  in  laminar  flows  as  well  as  a  jump  to  the 
turbulent  flow  regime,  due  to  the  employment  of  an  acoustic  field  has  direct  positive 
implications  on  heat  transfer  (Fraenkel  et  al.,  1998).  In  combustion  studies  acoustic  wave 
amplitude  has  been  correlated  to  heat  transfer  (Kwon  and  Lee,  1985).  temperature 
(Crocco,  1969),  emissions  (Keller  and  Hongo,  1990),  particle  dispersion  (Barrere  and 
Williams,  1969),  mixing  of  species  (Dubey  et  al.,  1996),  and  flame  structure  (McManus 
et  al.,  1990).    An  example  of  the  effect  of  acoustics  on  flame  structure  is  shown  in  Figure 
2.4  (Photographs  by  author).  These  photographs  of  propane  diffusion  flames  show 
marked  changes  in  the  flow  field  due  to  acoustic  oscillation.  Enhanced  mixing  is 
demonstrated  by  the  compacted  flame  zone  for  the  case  with  the  oscillating  flow  field.  A 
change  of  flow  field  direction  is  manifested  by  the  flame  zone  being  moved  downward 
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below  the  burner  arms.  Some  of  these  changes  seen  in  past  processes  can  be  expected 
with  implementation  of  an  acoustic  field  in  the  steam-reforming  process. 

i 


(A) 


(B) 


Figure  2.4:  A  15  kW  diffusion  propane  flame  without 
oscillating  flow  (A)  and  with  oscillating  flow  (B). 

Steam-Reforming  in  Combination  with  Acoustics 
Superimposing  an  oscillating  flow  in  the  axial  direction  of  the  reformer  (also 
referred  to  as  organ  pipe  oscillations)  with  the  steady  flow  component  is  proposed 
(Erickson  and  Roan,  2000).  High  acoustic  pressures  and  corresponding  local  particle 


22 
velocity  perturbations  are  achieved  by  establishing  resonance  in  the  reactor.  If  the 
acoustic  velocity  component  is  high  enough,  the  steady  flow  component  carries  the 
overall  flow  forward  and  the  oscillating  flow  portion  incrementally  steps  the  fluid 
forward  and  backward  across  portions  of  the  catalyst  bed.  This  allows  for  a  "mixing"  of 
the  boundary  layer  along  the  reactor  walls  and  past  the  catalyst  particles.  Under  the 
proper  conditions,  increased  heat  and  mass  transfer  rates  also  occur  as  a  result  of  the 
thinned  boundary  layers  and  acoustic-induced  mixing  in  the  reactor.    With  increased  heat 
and  mass  transfer,  temperature  gradients  are  also  reduced  throughout  the  acoustically 
enhanced  reactor.  The  projected  effects  of  acoustics  on  the  limiting  steps  of  heat  transfer, 
mass  transfer  and  chemical  kinetics  of  the  overall  process  are  discussed  as  potential 
theory-based  models  are  developed  below. 
Increased  Particle  Path  Length  due  to  Acoustics 

A  simplified  model  of  the  acoustics  shows  a  significant  increase  in  the  particle 
path  length,  resulting  in  a  decreased  effective  space  velocity,  without  changing  the  overall 
bulk  residence  time.  In  order  to  illustrate  this  effect,  consider  an  empty  open-ended  pipe 
(without  obstruction  from  the  catalyst).  Acoustic  displacement  amplitude  is  given  in 
Equation  2.9  (Halliday  and  Resnick,  1986),  where  Ym  is  the  maximum  acoustic 

ym  =  ^~  (2.9) 

displacement  amplitude  (m),  Pm  is  the  maximum  acoustic  pressure  (Pa),  f  is  frequency 
(Hz),  p  is  the  non-oscillating  density  of  the  gas  (kg/m3),  and  c  is  the  sonic  velocity  in  the 
gas  (m/sec).    The  equation  for  maximum  particle  velocity  due  to  acoustic  oscillation  (Um) 
(Kinsler  et  al.,  1982)  is  shown  in  Equation  2.10.  This  particle  velocity  is  subsequently 
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^m=—  (2.10) 

pc 

referred  to  as  acoustic  velocity.  In  order  to  derive  the  total  path  length  seen  by  a  particle 
traveling  through  the  reactor  (effective  particle  path  length),  one  must  consider  the  length 
of  the  reactor  (original  particle  path  length  with  no  acoustic  oscillation)  and  distance 
traveled  by  a  particle  due  to  acoustic  oscillation. 

A  simulation  of  the  particle  path  can  be  made  by  investigating  the  velocity 
components  during  the  residence  time.  The  slip  velocity,  or  velocity  relative  to  a 
stationary  catalyst  pellet  or  wall  is  the  sum  of  the  oscillating  velocity  component  Um 
(multiplied  by  a  spatial  factor  (Sf)  resulting  from  the  mode  shape  of  the  standing  wave) 
and  average  velocity  (Vs).  The  spatial  factor  is  assumed  to  follow  a  cycle  sine  wave 
pattern  within  the  reactor.  The  absolute  value  of  the  total  velocity  multiplied  by  the 
residence  time  yields  the  slip  displacement  (Xshp)  past  the  catalyst  or  wall.  If  the  time  is 
decreased  to  a  small  time  step  (At),  then  summed  over  the  residence  time  (tr)  the  slip 
displacement  can  be  represented  by  Equation  2.1 1 .  The  slip  displacement  can  be 
nondimensionalized  by  dividing  by  the  reactor  length.  Plotting  this  dimensionless  slip 
displacement  versus  the  velocity  ratio  (Um/Vs)  Equation  2.1 1  yields  the  total  theoretical 
slip  displacement  in  Figure  2.5 

Xs„P  =  i;;|(srUmsin27rft+Us)|At  (2.11) 

This  figure  shows  that  an  slope  asymptote  of  0.407  is  developed  at  large  velocity 
ratios.  This  limiting  slope  can  be  developed  through  investigation  of  the  displacement 
due  to  the  acoustic  field  within  the  reactor  during  the  residence  time.  As  before. 
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assuming  that  a  standing  wave  pattern  equivalent  to  a  single  cycle  sine  pattern  for 
velocity  is  developed  in  the  reactor,  a  spatial  average  displacement  amplitude  of  0.64  of 
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Figure  2.5:  Total  slip  displacement  and  slip  displacement  due  to  acoustics  for  increasing 

velocity  ratio. 

the  maximum  amplitude  will  occur.  The  equation  showing  the  displacement  of  the 
oscillating  portion  while  residing  in  the  reactor  is  shown  as  Equation  2.12  and  further 
simplified  in  Equation  2.13,  where  Xshp  is  the  slip  displacement  (m),  L  is  the  length  of  the 


X,.  n  =  (4)  x  — f—  x  0.64  x  —  f 
2Kfpoc  Vs 


LL  L 


Xslip  =  0.407  x 


(2.12) 


(2.13) 


reactor  (m),  and  Vs  is  the  steady  flow  component  of  the  reactor  (m/s).  One  can  see  that 
the  slip  displacement  is  independent  of  frequency  in  this  model.  This  results  from  the 
fact  that  while  an  increased  frequency  does  cause  more  oscillations  during  the  residence 
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time  the  amplitude  of  those  oscillations  is  decreased  proportionally  for  a  given  acoustic 
pressure. 

As  Figure  2.5  indicates,  the  model  predicts  large  increases  in  particle  path  length 
for  decreasing  average  velocity  or  increasing  particle  acoustic  velocity  by  higher  acoustic 
pressure  levels.  For  a  reasonable  example,  a  sound  pressure  level  of  165  dB  (re  20uPa)  in 
air  and  a  steady  flow  component  of  1  m/s  yields  a  particle  path  length  of  7  times  the 
particle  path  length  with  non-oscillating  flow.  If  employed  in  a  reforming  process  the 
effective  space  velocity  would  be  7  times  smaller  than  the  actual  space  velocity  or  the 
space  velocity  without  acoustics.  This  effect  is  diminished  with  the  longer  flow  path  and 
wave  attenuation  due  to  obstruction  from  catalyst  pellets.  Also  because  of  the  fact  that 
the  catalyst  is  not  continuous  throughout  the  catalyst  bed  these  results  must  be  scaled  by 
the  catalyst  packing  factor  and  a  factor  accounting  for  the  effective  area  of  the  catalyst 
bed  seen  by  the  acoustic  wave  (ie.  Exterior  of  the  catalyst  pellet).  However  this  model 
yielding  the  effective  space  velocity  gives  an  upper  boundary  of  the  expected  benefits  due 
to  an  acoustic  field  in  conjunction  with  methanol-steam  reformation  in  an  otherwise 
undisturbed  and  unobstructed  (without  catalyst  pellets)  flow  field. 
Increased  Heat  Transfer  Rate  due  to  Acoustics 

A  more  realistic  expectation  may  be  that  of  modeling  an  increase  in  capacity  due 
to  increased  heat  transfer  in  the  reformer.  As  presented  above,  heat  transfer  rate  is  seen  as 
a  fundamental  resistive  step  in  the  fuel  reformation  process.  Recent  papers  have 
addressed  this  issue  (Ohl  et  al.,  1996;  Dusterwald  et  al.,  1997;  Nagano  et  al..  2001)  and 
certain  methods  of  increasing  heat  transfer  rate  have  been  reported,  including  the  use  of 
fins  to  increase  the  heat  transfer  area  (Shiizaki  et  al.,  1999).  However,  increasing  heat 


26 


transfer  rates  with  acoustic  fields  for  other  applications  has  been  well  documented  over 
the  years  (Fraenkel  et  al.,  1998:  Erickson  et  al.,  1997:  Zinn,  1992:  Yavuzkurt  et  al.,  1991 ) 
Space  and  time  averaged  Nusselt  numbers  have  been  shown  by  Ha  and  Yavuzkurt 
(1993),  for  a  single  spherical  particle  in  a  pulsating  flow  with  a  steady  flow  component,  to 
be  represented  by  Equations  2.14  and  2.15.  In  these  equations,  Vr  is  velocity  ratio 

(2.15) 
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forV>\       -?-=  0.73839  V,       +0.16161 
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(previously  defined  as  acoustic  velocity  Um  divided  by  average  velocity  Vs),  Nup  is  the 
Nusselt  number  due  to  pulsating  flow,  and  Nus  is  the  steady  flow  Nusselt  number.  Nusselt 
number  is  a  direct  indicator  of  convective  heat  transfer  rate  as  shown  in  Equation  2.16, 
where  h  is  the  convective  heat  transfer  rate  coefficient,  k  is  the  thermal  conductivity  of 


hL 

Nu  =  — 
k 


(2.16) 


the  gas,  and  L  is  a  characteristic  length.  Figure  2.6  is  a  plotted  representation  of 
Equations  2.14  and  2.1 5.  While  interactions  between  particles  in  a  fixed  bed  reactor  are 


Figure  2.6:  Nusselt  number  ratio  versus  velocity  ratio  for  a  single  particle  in  an 

empty  duct. 
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clearly  much  more  complex  than  the  single  particle  example.  Figure  2.6  shows  the 
potential  for  increasing  convective  heat  transfer  rates  with  acoustic  fields.  For  the 
previous  example  of  a  SPL  of  165  dB  and  a  1  m/s  steady  or  average  velocity  component, 
the  convective  heat  transfer  rate  is  increased  by  a  factor  of  4  by  using  an  acoustic  wave. 
Fraenkel  et  al.  (1998)  have  found  similar  results  for  velocity  ratios  above  1  when  flow 
oscillations  are  used  in  grain  drying. 

The  increase  in  heat  transfer  rate,  by  definition,  allows  a  faster  response  to 
transients  including  load  following  and  start-up.  Mass  flow  capacity  can  also  been  seen 
to  increase  by  referring  to  a  simplified  energy  balance  within  the  catalyst  bed.  From  an 
energy  balance  in  the  reformer  (taking  into  account  convective  heat  transfer  from  the 
walls  to  the  gas  in  the  catalyst  bed),  allowable  reactant  mass  flow  rate  is  directly 
proportional  to  the  heat  transfer  rate  from  the  hot  wall  to  the  gas  as  shown  in  Equation 
2.17  where  m  is  mass  flow  rate,  h  is  the  convective  heat  transfer  rate  coefficient.  A  is  the 

•      hA(T    -  Tt  .    ) 

V      watt  local  uus  '  /■%    t  H\ 

m- —  (2.17) 

AH 

heat  transfer  area,  T  is  temperature,  and  AH  is  the  change  in  enthalpy  between  reactants 

and  products.  This  relation  applies  for  heat  transfer  limited  cases  that  are  not  limited  by 

diffusion  or  chemical  kinetic  rates. 

Increased  Overall  Reaction  Rate  due  to  Acoustics 

With  a  higher  heat  transfer  rate,  the  rate  limiting  steps  of  diffusion  and  chemical 

kinetics  become  more  important  to  the  reformation  process.  Although  mass  transfer  rate 

has  not  been  a  perceived  limitation  in  previous  work  (Asprey  et  al.,  1999:  Peppley  et  al.. 

1997).  smaller  reactors  with  high  llow  rates  will  deviate  from  the  plug  flow  reactor 
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assumption.  As  described  and  quantified  by  Ha  and  Yazavkurt  (1993),  acoustic  mixing 

increases  the  mass  flow  rate  at  the  surface  of  a  particle;  thus  in  a  reformer,  acoustics 

would  reduce  the  diffusion  time  of  a  mixture  to  reach  the  outside  of  a  catalyst.  Equations 

2.14  and  2.15  apply  to  the  Sherwood  number  ratio  as  well  as  to  the  Nusselt  number  ratio 

as  presented.  The  Sherwood  number  (Sh),  as  shown  in  Equation  2.18  is  analogous  to  the 

hmL 
Sh=-^-  (2.18) 

Nusselt  number  in  relation  to  convective  mass  transfer  rate.  In  Equation  2.1 8  hm  is  the 
convective  mass  transfer  rate  coefficient,  DAB  is  the  diffusivity  of  gas  A  into  the  gas  B, 
and  L  is  a  characteristic  length.    This  increase  in  convective  mass  transfer  should 
increase  the  diffusion  of  both  the  reactants  to  the  catalyst  particle  surface  and  the  products 
to  the  bulk  stream.  Internal  diffusion  or  mass  transfer  into  the  core  catalyst  area  is  not 
expected  to  be  significantly  affected  by  the  acoustic  waves  due  to  extremely  small 
dimensions  and  velocities,  thus  overall  diffusion  rate  enhancement  by  acoustics  may  be 
limited.  If  the  internal  diffusion  is  not  the  limiting  factor,  but  rather  the  process  is  limited 
by  external  diffusion,  acoustics  may  have  a  significant  effect. 

While  increased  mixing  from  acoustics  may  decrease  a  rate  limiting  diffusion  step 
to  a  catalyst  surface,  there  is  probability  that  the  chemical  kinetics  will  be  increased  as 
well.  Considering  Arrhenius  type  rate  equations,  reactions  will  proceed  more  quickly  for 
higher  temperatures.  With  the  reduction  of  local  temperature  gradients  due  to  acoustically 
induced  mixing  and  increased  heat  transfer,  the  reactor  controlled  by  local  temperature 
feedback  (as  discussed  by  Sterchi,  2001  and  Amphlett  et  al.,  1999)  would  have  higher 
temperatures  near  the  entrance  of  the  reactor  and  thus  higher  chemical  kinetic  rates  would 
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result.  Because  of  the  rate  limiting  step  of  methanol  decomposition  occurs  toward  the 
entrance  of  the  catalyst  bed.  minimization  of  the  temperature  gradient  at  this  point  is  vital. 
Additionally,  Tonkovich  et  al.  (1999)  have  found  that  for  the  water  gas  shift  reaction, 
intrinsic  reaction  times  can  be  on  the  order  of  milliseconds  rather  than  the  observed 
reaction  times  of  6-9  seconds.  They  exploit  the  fast  intrinsic  reaction  times  by  using 
micro-channel  reactors  where  fast  heat  and  mass  transfer  rates  govern.  In  a  similar 
approach,  acoustic  vibrations  should  increase  heat  and  mass  transfer  rates  at  the  exposed 
catalyst  surface  and  thus  may  decrease  the  required  contact  time  on  the  catalyst  surface. 

Although  it  is  obvious  that  not  all  of  the  area  in  a  catalyst  would  be  affected  by  the 
acoustic  wave  (due  to  internal  pores  etc.),  bounds  of  the  reaction  rate  can  be  established. 
The  lower  bound  would  be  the  case  where  the  reaction  rate  is  limited  only  by  internal 
diffusion  rate  and  thus  the  overall  reaction  rate  would  not  change  when  acoustics  are 
applied.  The  upper  bound  would  be  where  the  rate  is  wholly  limited  by  diffusion  rates  at 
the  surface  and  the  total  effective  area  (derived  from  the  effective  length)  would  be  used 
in  the  reaction  rate  models.  It  is  expected  that  with  catalyst  pellets,  internal  area  and 
diffusion  to  those  sites  should  dominate,  thus  deviating  from  the  upper  boundary 
predicted  by  use  of  the  effective  area.  With  screen  or  monolith  type  catalysts,  where 
acoustic  waves  would  have  more  of  an  effect  on  the  active  reaction  sites,  the  actual  rate 
may  follow  predicted  rate  using  the  effective  area  (derived  from  effective  length  in  Eq. 
2.1 1 )  more  closely. 
Extension  of  Catalyst  Life  by  Using  Acoustics 

In  some  aspects,  catalyst  life  is  also  expected  to  be  positively  affected  by  using  an 
acoustic  Held  in  the  reformation  process.  Acoustics  may  hinder  catalyst  degradation  due 
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to  poisoning,  fouling  and  sintering.  While  all  types  of  poisoning  may  not  be  affected  by 
the  acoustic  wave,  poisoning  by  oxygen  with  reduced  catalysts  could  be  mitigated  by  the 
oscillating  flow  under  the  proper  conditions.  With  an  oscillating  flow,  hydrogen 
produced  in  the  reformer  would  be  recirculated  to  the  upstream  catalyst  particles, 
effectively  re-reducing  any  oxidized  catalyst.  This  recirculating  characteristic  of  the  flow 
field  would  happen  at  velocity  ratios  (Vr)  above  1 .  Fouling  is  also  a  degradation 
mechanism  that  would  theoretically  be  minimized  by  employing  an  acoustic  field.  Much 
like  the  oscillations  of  a  sonic  cleaning  application,  clogging  at  the  surface  reaction  sites 
and  catalyst  pores  could  be  reduced  by  the  oscillating  action  of  the  passing  fluid. 
Reduction  of  certain  types  of  sintering  may  also  occur  as  spatial  temperature  gradients 
throughout  the  reactor  are  alleviated  and  temperature  fluctuations  are  minimized  by 
increased  control. 

There  also  may  be  a  mechanism  for  catalyst  degradation  due  to  acoustic  fields. 
Theoretically,  oscillating  pressure  fields  may  break  apart  the  catalyst  structure  or  drive 
particles  into  the  pores  of  the  catalyst  increasing  the  fouling  mechanism  or  creating  the 
effects  of  catalyst  sintering. 
Summary 

In  theory,  the  possible  benefits  which  arise  from  employing  acoustic  waves  in  the 
reforming  process  are  an  increase  in  capacity  for  a  given  catalyst  mass,  a  decrease  in 
transient  response  time,  an  increase  in  the  overall  efficiency  of  the  reformation  process, 
and  potential  extension  of  catalyst  life.  The  benefits  of  using  acoustics  in  related 
processes  including  combustion  reactors  have  been  shown  in  the  literature  as  increased 
convective  heat  transfer  rates,  increased  mass  transfer  rates,  and  increased  specie  mixing. 
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While  the  acoustic  field  may  not  significantly  affect  diffusion  and  chemical  kinetics  for 
all  catalysts,  especially  those  with  high  amounts  of  internal  area  not  affected  by  the 
acoustic  wave,  the  theoretical  investigation  points  toward  potential  for  enhancing 
reformation  processes  with  acoustic  waves. 

Contribution 
The  present  study  contributes  to  the  technical  field  in  that  it  is  both  a  theoretical 
development  as  presented  above  and  an  empirical  quantification  of  the  benefits  that  arise 
from  using  an  acoustic  field  in  conjunction  with  the  steam-reforming  of  a  hydrocarbon 
fuel.  As  shown  by  the  literature  review,  a  properly  applied  acoustic  field  can  benefit  other 
processes.  The  acoustic  field  in  the  reformation  process  results  with  similar  enhancement 
and  is  expected  to  allow  smaller  size  and  lower  weight  of  the  reforming  reactor,  which  in 
turn  would  yield  quicker  start-up  times  and  faster  dynamic  response.  The  lifetime  of  the 
catalyst  is  expected  to  be  extended  by  utilizing  acoustic  waves.  The  present  investigation 
quantifies  these  benefits  and  may  contribute  towards  a  practical  method  of  reforming  a 
liquid  fuel  with  high  efficiency  for  fuel  cell  vehicles. 


CHAPTER  3 
EXPERIMENTAL  APPROACH  AND  FACILITY 

This  chapter  discusses  the  experimental  approach  and  facility.  A  description  of 
the  methanol  steam  reforming  facility  is  presented.  Various  subassemblies  of  the  facility, 
as  well  as  instrument  and  control  interfaces  are  described.  This  includes  discussion  of  the 
acoustic  field  generator  and  schemes  for  controlling  the  temperature  and  acoustic  field. 
The  instrumentation  is  reviewed  along  with  the  methods  used  for  collecting  and  analyzing 
data.  Modeling  of  the  reformation  process  is  discussed,  which  discussion  includes  the  set 
up  of  a  factorial  experimental  design. 

Experimental  Approach 

The  experimental  approach  was  designed  to  quantify  the  effect  of  the  acoustic 
field  on  the  reformation  process.  This  included  adaptation  of  a  reformer  facility  to  accept 
acoustic  waves,  development  of  procedures  for  operation  of  the  facility,  establishment  of 
temperature  and  acoustic  control  in  the  catalyst  bed,  the  acoustical  analysis  of  that 
facility,  and  the  design  of  an  experiment  from  which,  a  statistically  valid  empirical  model 
could  be  developed.  Further  data  were  taken  to  show  curvature  of  the  model,  verify  the 
model  results,  demonstrate  the  effect  of  acoustics  on  catalyst  degradation,  power  demand, 
as  well  as  on  the  reactor  temperature  profile  and  controllability.    Analysis  of  the  data 
enabled  construction  of  a  hydrogen  power  output  per  acoustic  power  input  ratio  which 
quantifies  the  energy  benefit/cost  on  the  reformation  process  from  the  acoustic  field. 
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Experimental  Facility 
A  reforming  facility  was  previously  constructed  at  the  University  of  Florida  Fuel 
Cell  Research  and  Training  Laboratory.  This  reforming  facility,  shown  in  Figure  3.1 
(photograph  by  J.  Patrick  Sterchi)  has  been  used  in  the  past  to  study  the  effects  of 


Figure  3.1 :  Reforming  facility  at  the  University  of  Florida  Fuel  Cell  Research  and 

Training  Laboratory. 

impurities  in  methanol  on  the  steam-reforming  process  (Sterchi.  2001 ).  The  reforming 

facility  hardware  consists  of  pumps,  vaporizers,  superheaters,  two  electrically  heated 

reforming  reactors,  a  condensing  unit,  a  power  relay  box.  and  thermal  instrumentation  in 

connection  with  computer  control  for  the  power  relay  box  and  pumps.  As  the  facility  was 
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originally  designed  to  simultaneously  operate  two  reforming  reactors,  the  present  study 
utilized  only  one  side  of  the  total  facility.  Figure  3.2  is  a  simplified  schematic  of  the 
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Figure  3.2:  Simplified  schematic  of  the  methanol  steam  reforming  processor, 
system  used.  In  this  system  methanol :water  premix  is  pumped  through  an  electrically 
heated  vaporizer.  The  vaporized  gas  is  routed  through  an  electrically  heated  liquid  trap  or 
secondary  vaporizer.  This  liquid  trap  insures  complete  vaporization  of  the  premix 
entering  the  reformer  at  high  flow  rates.  The  gaseous  premix  then  may  be  routed  to  the 
exhaust  or  directed  to  the  superheater  by  manipulation  of  the  valves  shown.    If  routed  to 
the  electrically  heated  superheater  the  premix  follows  the  path  shown  to  the  catalyst  bed. 
An  oscillating  pressure  is  directed  through  the  catalyst  bed  from  a  nitrogen-purged  and 
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water-cooled  acoustic  driver  attached  at  the  lower  end  of  the  catalyst  bed.  After  traveling 
through  the  electrically  heated  catalyst  bed  the  reacted  mixture  can  be  routed  through  the 
condensing  unit  or  to  the  exhaust  by  manipulation  of  the  valves  shown.  After  passing 
through  the  condensing  unit  the  dry  gas  may  be  sampled  through  a  gas  sampling  port  as  it 
travels  to  the  exhaust.  A  carbon  dioxide  purge  unit  attaches  to  the  superheater  assembly 
to  keep  the  superheater  and  catalyst  bed  free  from  oxygen  during  non-operating  periods. 

The  overall  system  was  divided  into  separate  subassemblies.  The  methanol: water 
premix  reservoir,  pump,  and  flow  metering  devices  were  grouped  into  the  pumping 
subassembly.  The  vaporizer,  liquid  trap  (or  secondary  vaporizer)  and  the  superheater 
comprise  the  vaporizer  subassembly.  The  reactor  subassembly  consists  of  the  catalyst  bed 
and  acoustic  components.  These  acoustic  components  are  the  acoustic  driver,  a  nitrogen 
purge,  a  water  cooling  system,  a  pressure  transducer,  and  adapters.  The  condensing  unit 
and  gas  sample  port  make  up  the  condenser  subassembly.  Computer  control  and  logic  for 
both  temperature  and  acoustic  field  regulation  is  explained  in  the  control  subassembly. 
Each  subassembly  is  more  fully  described  as  a  separate  section  below. 
Pumping  Subassembly 

The  pumping  subassembly  consists  of  a  reservoir  placed  on  a  scale,  a  peristaltic 
pump  and  associated  tubing  as  shown  in  Figure  3.3.  The  scale  and  a  timing  device  was 
used  to  calculate  the  actual  mass  flow  rate  of  premix  fed  into  the  reactor.    The  scale  used 
was  an  Acculab  VI-4kg  that  had  a  resolution  of  1  gram.    The  peristaltic  pump  was  a 
Masterflex  L/S  Model  7519-20  cartridge  pump.  Rough  control  of  premix  mass  flow  rate 
was  enabled  by  changing  the  volumetric  flow  setting  on  this  peristaltic  pump.  The 
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volumetric  flow  setting  could  be  set  at  increments  of  0.1  ml/min.  The  actual  mass  flow 
rate  was  found  by  measuring  the  mass  change  of  the  reservoir  per  unit  time. 
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Figure  3.3:  Artistic  representation  of  the  pumping 
subassembly. 

Vaporizer  Subassembly 

The  vaporizer  subassembly  consisted  of  a  vaporizer  and  a  liquid  trap  or  secondary 
vaporizer  as  shown  in  Figure  3.4,  and  the  superheater  shown  in  Figure  3.5.  The  vaporizer 
was  a  1 .27  cm  (0.5  in.)  O.  D.  stainless  steel  tube  cut  to  85  cm  (33.5  in.)  in  length.  The 
length  of  the  tube  was  wrapped  in  a  244  cm  (96  in.)  650  Watt  heat  tape.  The  liquid  trap 
or  secondary  vaporizer  was  also  made  of  1.27  cm  (0.5  in.)  O.  D.  stainless  steel  tubing 
that  was  cut  to  64  cm  (25.2  in).  This  tube  was  wrapped  in  a  1 83  cm  (72  in.)  length  435 
Watt  heat  tape.  The  temperatures  at  the  pipe  surfaces  and  the  exits  were  monitored  by 
stainless  steel  sheathed  k-type  thermocouples.  The  superheater  assembly  as  shown  in 
Figure  3.5  is  a  2.68  cm  O.D.  (nominal  3/4"  schedule  40)  stainless  steel  pipe  cut  to  a  38 
cm  (15  in.)  length.  Receptacles  for  fittings  (7/8  in. -14)  were  machined  into  both  ends  of 
the  pipe.  The  pipe  was  wrapped  in  a  366  cm  (144  in.)  430  Watt  ceramic  beaded  heater 
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Figure  3.4:  Artistic  representation  of  the  vaporizer  and  liquid  trap, 
and  outer  surface  temperature  was  measured  near  the  bottom  of  the  assembly  by  a 
sheathed  k-type  thermocouple.    The  exit  of  the  pipe  was  made  by  allowing  a  15.25  cm 
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Figure  3.5:  Artistic  representation  of  the  superheater 

(6  in.)  long  0.635  em  (0.25  in.)  O.D.  stainless  steel  tube  to  be  placed  at  the  bottom  of  the 
assembly  recessing  through  the  lower  fitting  into  the  superheater  cavity.    This  recessed 
outlet  insured  that  no  liquid  phase  fluid  accumulating  at  the  bottom  of  the  superheater 
would  be  allowed  into  the  catalyst  bed.  The  exit  temperature  was  monitored  at  the  exit  of 
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the  recessed  outlet  pipe  by  a  stainless  steel  sheathed  k-type  thermocouple  inside  the 
fitting  directly  above  the  catalyst  bed. 
Reactor  Subassembly 

The  reactor  subassembly  consisted  of  the  catalyst  bed.  and  acoustic  components. 
These  two  sections  are  discussed  individually  below. 
Catalyst  bed 

As  shown  in  Figure  3.6,  The  catalyst  bed  was  housed  in  a  2.68  cm  O.D.  (nominal 
3/4  in.  schedule  40)  stainless  steel  pipe  cut  to  a  63.5  cm  (25  in.)  length.  Receptacles  for 
fittings  were  made  at  both  ends  of  the  pipe.  The  bed  was  filled  with  a  pelletized  low- 
temperature  water-gas  shift  catalyst  manufactured  by  ICI  Katalco  (ICI  catalyst  51-3).  The 
catalyst  was  designed  for  a  temperature  range  of  200-300 °C.  This  catalyst  is  a 
proprietary  combination  of  copper  oxide,  zinc  oxide,  and  aluminum  oxide  and  had  to  be 
reduced  before  use.  The  reduction  process  for  this  catalyst  takes  place  over  40  hours  at 
1 80 °C  with  a  2.3%  hydrogen,  97.7%  nitrogen  gas  mixture,  at  a  gas  hourly  space  velocity 
of  300  hr"1.  The  40  hour  period  is  followed  by  an  hour  at  210°C  with  a  6.5%  volume 
fraction  of  hydrogen.  The  reduction  finishes  with  a  final  hour  at  210°C  with  a  12  % 
volume  fraction  of  hydrogen.  This  reduction  process  is  explained  by  Sterchi  (2001 ).  The 
pelletized  catalyst  was  dimensionally  0.52  cm  in  length  and  0.54  cm  in  diameter  and  had 
a  bulk  density  of  1 .2-1 .3  kg/liter.  In  order  to  avoid  blockage  of  the  exit,  the  catalyst  bed 
was  supported  by  a  stainless  steel  screen  fixed  at  the  top  of  the  adapter  for  the  acoustic 
generator.    The  catalyst  bed  was  divided  into  6  zones.  At  the  end  of  each  zone  and  before 
the  first  zone  a  0.93  cm  diameter  (1/8  in.  NPT)  port  was  machined  into  one  or  both  sides 
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Figure  3.6:  Artistic  representation  of  the  catalyst  bed  housing, 
of  the  stainless  steel  pipe.  Fittings  could  be  placed  into  the  machined  holes  to  extract 

temperature  or  pressure  information  at  each  respective  location  in  the  bed. 

The  fittings  were  sealed  by  a  high  temperature  silicone  sealant.  Above  the  first 

zone,  two  ports  are  used  for  monitoring  dynamic  pressure  and  monitoring  the  exit 
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temperature  of  the  superheater.  Each  zone  was  individually  heated  by  a  250  Watt  band 
heater.  Control  for  these  heaters  was  determined  by  the  temperature  at  the  respective 
zone  exit  as  explained  in  the  control  subassembly  below.  Temperature  was  monitored 
using  k-type  stainless  steel  sheathed  thermocouples.  Each  fitting  was  machined  to  be 
flush  with  the  interior  surface  of  the  pipe  at  the  zone  exits.  The  control  thermocouples 
were  aligned  as  to  be  near  the  interior  wall  location  in  the  fitting.  These  temperatures 
represented  the  interior  wall  temperature  at  the  specified  axial  location.    At  the  bottom  of 
the  pipe  an  adapter  was  threaded  into  the  machined  receptacle  to  allow  an  acoustic  driver 
to  modulate  the  catalyst  bed  pressure. 
Acoustic  components 

Transmission  of  acoustic  pressure  and  associated  velocity  fluctuations  are  a 
product  of  geometric  conditions,  and  also  the  fluid  properties.  Additionally,  transmission 
is  typically  inhibited  by  obstruction  in  the  acoustic  path.  In  a  packed  catalyst  bed,  the 
catalyst  itself  partially  obstructs  the  propagation  of  acoustic  waves.  Notwithstanding  this 
difficulty,  standing  acoustic  waves  can  build  up  within  the  bed  as  a  result  of  wave 
reflection  at  the  ends  of  the  reactor. 

An  acoustic  field  was  generated  by  pulsating  an  acoustic  driver  attached  to  an 
adapter  at  the  bottom  of  the  reactor  as  shown  in  Figure  3.7.  The  acoustic  driver  used  was 
a  Klipsch  Pro  1 133  Horn  Driver  typically  used  in  amphitheater  applications.  The  method 
of  using  loudspeakers  to  modulate  pressure  in  a  chemical  reaction  has  been  proven  by 
several  researchers  including  McManus  et  al.  (1990),  Poinsot  et  al.  (1989).  Heckl  (1988), 
Lang  et  al.  (1987).  Zikikiout  et  al.  (1986),  and  many  others.  The  horn  driver  was  attached 
to  mounting  devices  and  directed  into  the  catalyst  bed  as  shown.  As  the  horn  driver  was 
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Figure  3.7:  Artistic  representation  of  the  acoustic  driver,  driver  housing,  mounting  plate 

and  fittings. 

not  designed  to  be  used  in  relatively  high  temperature  chemical  applications,  a  water 

cooling  jacket  was  designed  around  a  fitting  to  eliminate  conduction  from  the  catalyst  bed 

walls.  This  water-cooled  section  was  also  a  heat  sink  for  the  reactor  as  it  pulled  heat 

away  from  the  bottom  zone.  A  nitrogen  purge  was  directed  to  both  the  front  of  the 

diaphragm  and  the  back  of  the  driver  housing.  The  nitrogen  at  the  front  of  the  diaphragm 

was  used  to  inhibit  hot  reacted  species  from  flowing  into  the  driver.  The  nitrogen  purge 

at  the  back  of  the  housing  was  a  safety  precaution  as  there  was  possibility  of  electrical 

ignition  in  the  voice  coils  in  the  presence  of  hydrogen.  The  entire  housing  and  fittings 

were  tilted  upward  at  a  10°  angle  to  allow  any  condensed  species  forming  on  the  water 

cooled  section  to  flow  down  toward  the  exit  of  the  reactor  rather  than  condensing  and 


43 
pooling  in  the  driver  or  fittings.  The  driver  housing  and  fittings  were  sealed  with  the  same 
silicone  sealant  as  the  port  fittings  in  the  high  temperature  zone. 

The  acoustic  wave  was  monitored  by  a  Kistler  piezoelectric  A206  pressure 
transducer.  This  piezoelectric  transducer  used  a  Kistler  piezotron  coupler  and  was  factor}' 
calibrated  at  15.3  mV/kPa.  As  previously  shown  in  Figure  3.6,  the  pressure  transducer 
was  mounted  at  the  port  above  the  first  zone  in  the  catalyst  bed  for  most  of  the 
experiments  performed.  This  position  was  selected  so  that  the  acoustic  wave  would  have 
to  travel  through  the  entire  length  of  the  catalyst  bed.  The  mounting  took  place  through  a 
1.0  cm  O.D.  (1/8  in.  schedule  40)  stainless  steel  pipe  cut  to  a  12.7  cm  (5  in.)  length  and 
an  adapter  to  connect  the  pressure  transducer.  A  water  cooling  jacket  was  constructed  for 
the  mounting  section  but  after  monitoring  temperatures  at  the  pressure  transducer  it  was 
determined  that  water  cooling  was  not  necessary.  Preliminary  experimentation  found  that 
water  cooling  at  the  transducer  condensed  liquid  species  directly  above  the  catalyst  bed. 
Experiments  performed  for  catalyst  degradation  required  the  elimination  of  this 
condensation  and  possible  intrusion  thus  the  pressure  transducer  was  removed  for 
experiments  dealing  with  degradation.  When  the  pressure  transducer  was  absent  in  the 
degradation  experiments,  acoustic  driver  feedback  came  from  the  voltage  input  to  the 
acoustic  driver. 
Condensing  Unit  Subassembly 

A  condensing  unit  was  also  used  that  cooled  the  exit  gases  from  approximately 
250°  C  to  0°  C.  As  shown  in  Figure  3.8,  This  unit  consisted  of  a  coiled  0.635  cm  (0.25 
in.)  O.  D.  copper  tubing  through  which,  the  exit  gases  flowed.  The  exterior  of  this  copper 
tubing  was  Hooded  with  water  from  an  ice  bath.  Reacted  species  would  condense  in  the 
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Figure  3.8:  Artistic  representation  of  the  condenser  subassembly, 
tubing  and  follow  the  flow  path  to  a  liquid  trap  where  the  liquid  species  would  remain 

and  the  dry  product  gases  would  be  directed  to  the  exhaust.  The  above  condensing 

components  were  housed  in  a  refrigeration  unit  that  decreased  the  heat  losses  of  the  ice 

bath.  Temperatures  were  monitored  at  the  dry  gas  exit  and  at  the  interior  of  the 

refrigeration  unit  with  sheathed  k-type  thermocouples.  A  dry  gas  sampling  port  was 

available  as  the  dry  gas  products  were  directed  to  the  exhaust. 

Instrumentation  and  Control  Subassembly 

Control  of  the  temperatures  and  acoustic  pressure  took  place  in  the  form  of  an 

active  closed  loop  and  open  loop  schemes,  respectively.  The  input  and  output  signals  for 

these  two  control  schemes  are  found  in  Figure  3.9.  The  instrumentation  for  these  systems 

were  the  output  signals  coming  from  the  reactor.  A  description  of  the  control  schemes  is 

given  in  the  following  subsections. 
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Figure  3.9:  Input  and  output  signals  for  temperature  and  acoustical  control. 


Temperature  control 

The  heating  elements  discussed  previously  were  computer  controlled  by  feedback, 
from  the  thermocouple  outputs.    The  14  thermocouple  signals  were  processed  using  a 
Iotech  DBK-52  Thermocouple  module  that  was  linked  to  a  data  acquisition  and  control 
program  that  was  built  by  the  author  in  LabVIEW.  This  temperature  control  program  can 
be  found  in  Appendix  A.  Control  of  the  9  heating  elements  took  place  through  an  Iotech 
DBK-24  Digital  Output  module  linked  to  the  same  program.  Depending  on  the 
individual  location  and  temperature  history  sensed  by  the  respective  thermocouples,  the 
output  module  switched  individual  relays  that  in  turn  switched  120  V  to  the  heating 
elements  through  a  solid  state  relay  as  shown  in  the  control  circuit  displayed  in  figure 
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3.10.  The  relays  were  contained  in  a  power  distribution  box  exhibited  in  Figure  3.11.  As 
shown  by  the  control  logic  diagram  found  in  Figure  3.12,  control  occurred  as  the 
thermocouple  signals  were  compared  to  a  setpoint  temperature.  The  slope  of  the 
temperature  signals  was  also  monitored  and  compared  to  predetermined  set  points  as  was 
the  duration  of  the  on  position  for  each  heating  element.  The  ability  to  override  the 
control  was  given  to  the  user  but  this  option  was  not  used  for  any  of  the  experiments 
presented.  The  setpoint  matrix  used  for  the  automatic  control  of  the  temperatures  is 
found  in  Table  3.1.  The  temperature  control  scheme  with  the  matrix  in  Table  3.1  was 
found  to  regulate  temperature  more  precisely  than  other  control  schemes  previously 
utilized  with  this  steam  reforming  reactor.  While  previous  studies  did  not  require  precise 
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Figure  3.10:  Control  circuit  for  heating  elements. 


47 


# 


Figure  3.11:  Power  distribution  box. 
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Table  3.1 :  Control  setpoint  matrix. 

Setpoint 
°C 

DT/dt 
(°C/sec) 

Maximum  Duration 
(%  time) 

Vaporizer  Surface 

420 

0.3 

100 

Vaporizer  Exit 

290 

20 

100 

Trap  Surface 

480 

0.3 

100 

Trap  Exit 

250 

20 

100 

Superheater  Surface 

305 

0.2 

100 

Superheater  Exit 

250 

0.1 

100 

Zone  1  Exit 

250 

0 

100 

Zone  2  Exit 

250 

0 

100 

Zone  3  Exit 

250 

0 

100 

Zone  4  Exit 

250 

0 

100 

Zone  5  Exit 

250 

0 

100 

Zone  6  Exit 

250 

0 

100 

control  of  the  instantaneous  temperature,  it  was  expedient  in  the  present  study  to 
eliminate  possible  variations  to  clearly  identify  factors  affecting  the  reforming  process. 
Figure  3.13  shows  15  minute  time  scans  of  the  Zone  1  exit  temperatures  with  a  control 
scheme  previously  developed  by  Sterchi  (2001 )  and  the  control  scheme  used  in  the 
present  study.  Both  temperature  plots  were  acquired  under  similar  operating  conditions 
with  a  control  setpoint  of  250°C  with  no  acoustic  wave,  and  a  premix  flow  rate  of  4.78 
g/min.  The  major  differences  between  these  two  control  schemes  are  the  control  matrix, 
the  location  of  the  feedback  thermocouples,  and  the  use  of  the  liquid  trap  (or  secondary 
vaporizer)  for  the  present  study.  Previous  operation  utilized  feedback  thermocouples  at 
the  radial  centerline  of  the  catalyst  bed  while  the  control  scheme  used  in  the  present  study 
obtains  feedback  from  thermocouples  at  the  interior  wall  positions.  This  position  near  the 
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Figure  3.13:  Time  plot  of  catalyst  bed  temperatures  at  zonel  exit  utilizing  two 
control  schemes  with  a  250  °C  setpoint. 

wall  allows  a  faster  response  time  to  the  heat  input  and  thus  the  overshoot  in  temperature 

is  reduced.  The  fluctuations  in  temperature  over  the  period  shown  in  Figure  3.13  for  the 

present  study  are  representative  of  temperature  measurements  over  an  entire  experiment 

period  and  throughout  the  catalyst  bed.  These  temperatures  and  fluctuations  are  also 

characteristic  for  all  the  experiments  performed  in  this  study.  The  average  and  standard 

deviation  for  all  zone  exit  temperatures  under  the  conditions  described  above  are  shown 

in  Table  3.2. 

Table  3.2:  Average  and  standard  deviation  for  all  catalyst  bed  temperatures  utilizing  a 

250°C  setpoint. 


Average  Temperature 

Standard  Deviation 

Zone  1  Exit 

251.4 

1.6 

Zone  2  Exit 

251.2 

1.5 

Zone  3  Exit 

250.9 

0.9 

Zone  4  Exit 

251.2 

1.3 

Zone  5  Exit 

254.1 

1.0 

Zone  6  Exit 

251.X 

3.4 
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Acoustic  control 

An  acoustic  field  was  generated  with  an  open  loop  control  scheme.  As  previously 
shown  by  the  input  and  output  signals  in  Figure  3.8,  a  monotone  sine  wave  or  swept  sine 
wave  reference  signal  was  induced  with  a  Wavetek  Sweep  function  generator.    This 
reference  wave  then  was  directed  to  the  channel  input  of  a  QSC  Mx-700  amplifier.    The 
amplifier  was  in  turn  connected  to  the  Klipsch  Pro  1 133  Horn  Driver  through  a  pressure 
housing.  The  amplified  AC  voltage  seen  by  the  driver  was  monitored  by  connecting  a 
handheld  Fluke  multimeter  in  parallel  with  the  driver.  The  acoustic  output  was  probed  by 
observing  the  conditioned  output  voltage  from  the  Kistler  piezotron  coupler  that  was 
connected  to  the  piezoelectric  pressure  transducer  as  described  previously.  Both  the  input 
(from  the  function  generator)  and  the  output  signal  (from  the  piezotron  signal 
conditioner)  were  simultaneously  monitored  by  an  analog  oscilloscope  and  a  digital  data 
acquisition  unit  programmed  in  LabVIEW  by  the  author.  These  digital  data  acquisition 
programs  can  be  found  in  Appendix  A.  These  programs  included  a  spectral  analysis 
program  for  measuring  the  acoustic  properties  of  the  reactor  via  transfer  functions  and 
anoscilloscope  program  for  measuring  and  recording  the  sound  pressure  level  data.  The 
frequency  of  oscillation  was  held  constant  throughout  the  experiment.  The  resonant 
frequencies  and  acoustic  field  properties  within  the  catalyst  bed  are  presented  and 
discussed  in  Chapter  4.  As  shown  in  Figure  3.14,  which  is  a  typical  time  trace  of  the 
acoustic  pressure,  there  were  both  fluctuations  in  the  acoustic  pressure  amplitude  levels 
and  60  Hz  noise  present  in  the  transducer  output  signal.  Because  the  noise  fluctuations 
increased  the  observed  maximum  amplitude,  the  average  noise  level  amplitude  found 
after  an  experiment  run  was  removed  from  the  average  maximum  found  during  the 
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Figure  3.14:  Typical  time  trace  of  pressure  in  the  catalyst  bed  while  generating  hydrogen 

with  an  acoustic  input  of  650  Hz  . 

experiment  run.  This  post  processing  yielded  a  more  accurate  representation  of  maximum 

acoustic  pressure  and  the  corresponding  sound  pressure  level  found  inside  the  reactor. 

The  sound  pressure  level  (SPL)  in  dB  was  found  by  applying  Equation  3.1  where  Pmis  is 

the  root  mean  square  of  the  pressure  wave  amplitude  and  Pref  is  a  standard  reference  of 

20uPa  used  throughout  the  results. 


SPL=  20  log 


'P     N 

nns 


vPrc,-; 


(3.1) 


Analysis  of  Data 
There  were  multiple  steps  and  sub-processes  required  for  analysis  of  the  collected 
data.  The  sub-processes  required  for  analysis  were  use  of  gas  chromatography, 
calculations  of  conversion,  and  calculations  of  space  velocity. 
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Gas  Chromatography 

An  SRI  model  8610-C  Gas  Chromatograph  (GC)  was  used  to  analyze  dry 
products  of  the  reactor.    The  GC  was  computer  controlled  using  Peak  Simple  software  as 
provided  by  the  manufacturer.  The  GC  was  equipped  with  a  91 .44  cm  (36  in.)  molecular 
sieve-packed  column  and  a  182.88  cm  (72  in.)  silica  gel-packed  column.  Tedlar  bags  (1.5 
liter)  were  used  to  collect  dry  gas  samples.  After  the  tedlar  sampling  bag  was  connected 
to  the  sampling  valve  on  the  gas  chromatograph,  a  40-50  ml  gas  sample  was  drawn 
through  the  sampling  valve  with  a  syringe.  The  automated  gas  sampling  valve  was  kept 
heated  to  100°C  and  was  used  to  inject  a  1  ml  sample  into  the  columns  upon  command 
from  the  computer  user  interface.  A  Bourdon  type  pressure  gauge  was  used  in  connection 
with  the  gas  sampling  valve  to  ensure  that  the  sample  in  the  valve  was  near  ambient 
pressure  before  the  sample  was  injected.  The  flow  paths  for  the  GC  sampling  valve, 
columns  and  detectors  are  shown  in  Figure  3.15.  As  the  gas  sampling  valve  would  rotate 
from  inject  back  to  load  position  (at  9.0  min).  the  flow  in  the  columns  would  reverse. 
This  allowed  separation  of  H2  02  N2  CH4  and  CO  specie  passing  through  the  mol  sieve 
and  separation  of  C02  and  higher  hydrocarbons  trapped  in  either  the  molecular  sieve- 
packed  column  or  silica  gel-packed  column  after  the  flow  reversal  (SRI,  1997). 
Separation  occurred  and  the  silica  gel-packed  column  was  heated  (molecular  sieve- 
packed  column  was  kept  at  100°C  for  all  experiments)  as  shown  in  Figure  3.16.  The 
detectors  used  were  a  thermal  conductivity  detector  (TCD)  and  a  tlame  ionization 
detector  (FID)  in  conjunction  with  a  methanizer.  The  voltage  traces  from  these  detectors 
and  a  calibration  file  were  used  by  the  software  to  calculate  the  volume  percent  of  species 
in  the  gas  sample.  The  TCD  was  used  to  detect  H2  02  and  N2.  as  well  as  high  levels 
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Figure  3.15:  Flow  paths  in  the  gas  chromatograph. 
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igure  3. 16:  Siliea  gel-packed  column  temperature  and  approximate  specie  elution  times 
for  the  species  expected  with  the  gas  chromatograph. 
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(above  1%)  of  CH4,  CO,  and  C02.  The  FID  was  utilized  to  detect  low  levels  of 
carbonaceous  compounds  including  CH4,  CO,  C02  and  non-methane  hydrocarbons 
(NMHC).  All  of  the  above  specie  (except  hydrogen)  were  detected  using  a  ultra  high 
purity  grade  helium  carrier  gas  and  a  high  TCD  current  setting.  Detection  of  H2  levels 
with  the  GC  required  use  of  an  ultra  high  purity  grade  nitrogen  carrier  gas  and  a  low  TCD 
current  setting. 

Calibration  of  the  GC  took  place  as  standardized  gas  samples  were  tested  in  the 
GC.  The  calibration  levels  for  each  gas  is  shown  in  Table  3.3.  With  the  exception  of  the 
higher  levels  of  hydrogen,  these  calibration  mixtures  were  acquired  from  calibrated 
sample  gas  cylinders. 

Table  3.3:  Gas  chromatograph  calibration  levels. 


Specie 

H2 

02 

N2 

CH4 

CO 

C02 

NMHC 

Calibration 
levels 

1%,4%,  10%, 

-65%,  -70%, 

-75%,  -80%, 

100% 

5%, 
21% 

5%, 
20%, 
79% 

20 
PPM, 

4% 

100  PPM, 
5% 

10%, 
20%, 
25% 

Ethylene 
20  PPM 

Rough  calibration  of  the  higher  hydrogen  volume  percentages  was  made  by  filling 
a  tedlar  bag  with  metered  hydrogen  and  nitrogen.  As  there  was  some  fluctuation  in  the 
metered  values  and  these  points  were  not  independently  verified,  this  calibration  was 
considered  approximate.  This  rough  calibration  did  however  verify  that  the  balance  of 
products  not  acquired  using  the  helium  carrier  was  hydrogen.  This  is  shown  by  an 
example  of  the  volume  percentages  found  for  the  dry  product  gases  during  preliminary 
data  at  various  operating  conditions  in  Table  3.4  below.  Table  3.4  also  demonstrates  the 
variability  due  to  the  gas  chromatograph  sampling  as  Experiments  1  and  2  were  gathered 
under  the  same  operating  conditions.  Because  of  this  demonstrated  variability  in  overall 
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species  concentration,  the  overall  average  concentration  of  hydrogen  was  used  to 
calculate  hydrogen  output  power  in  the  analysis.  Thus  the  hydrogen  output  power  was  a 
function  only  of  the  dry  product  gas  flow  rate  rather  than  being  affected  by  the  variable 
concentrations  during  a  particular  experiment. 

Table  3.4:  Volume  percent  species  of  the  dry  product  outlet  gas. 


Specie 

H: 

0: 

N2 

CH4 

CO 

co2 

NMHC 

Total 

Experiment  1 

68.76 

ND 

6.52 

ND 

0.54 

22.83 

ND 

98.65 

Experiment  2 

72.65 

ND 

4.04 

ND 

0.55 

23.92 

ND 

101.16 

Experiment  3 

71.19 

ND 

5.21 

ND 

0.50 

23.11 

ND 

100.01 

ND=Not  Detected 


Conversion 


Conversion  of  the  methanol  is  defined  as  shown  in  Equation  3.2  where  X  is  the 
conversion  of  the  methanol  fuel,  and  methanol  is  determined  on  a  mass  basis.  The 
calculated  conversion  is  then  independent  of  the  sample  time  period.  The  time  period  and 


y 

■^  CHjOH 


CHpHinpul  -  CHpHQulput 


CHpti 


(3.2) 


Input 


corresponding  amounts  of  condensate  did,  however,  affect  the  variability  associated  with 
the  conversion.  Shorter  time  periods  yielded  less  condensate  and  thus  these  small 
amounts  of  condensate  were  more  prone  to  contamination  from  water  droplets  that 
formed  on  the  liquid  trap  when  exposed  to  the  atmosphere  prior  to  analysis  or  water  that 
would  leak  into  the  liquid  trap  upon  opening.  It  was  determined  from  comparing  small 
amounts  of  condensate  from  shorter  experiment  times  to  those  of  longer  experiments 
under  the  same  conditions  that  a  minimum  of  40  grams  of  condensate  should  be  collected 
to  insure  negligible  biasing  of  conversion  data  by  sample  collection  time.    The  conversion 
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equation  shown  above  in  Equation  3.2  required  that  the  methanol  mass  fraction  in  the 
condensate  be  known.  The  amount  of  methanol  in  a  condensate  sample  was  determined 
by  utilizing  the  condensate  density.  The  condensate  density  was  found  by  way  of  a  Anton 
Paar  density  meter  that  had  a  0.0001  g/cm3  resolution.  The  density  (p)  was  input  to  an 
emperical  formula  developed  by  Sterchi  (2001)  that  yields  methanol  mass  fraction 
(mfCH30H)  of  the  condensate.    This  empirical  relation  is  found  in  Equation  3.3  and  was 
verified  for  this  study.  From  the  methanol  mass  fraction  of  the  condensate,  the  mass  of 
unused  methanol  (output)  is  found,  and  the  conversion  can  then  be  determined. 

mfcH3oH  =  3859.9p2  -  823  Up  +  4370.0  (3.3) 

Space  Velocity 

Space  velocity  for  this  investigation  was  calculated  from  the  mass  flow  rate  of 
methanol: water  premix.  Flow  rate  data  was  gathered  as  masses  of  both  premix  consumed 
and  condensate  accumulated  could  be  known  over  a  unit  of  time.  Space  velocity  is  found 
by  measuring  the  volume  of  the  reactor  (measured  catalyst  bed  length  was  used  for 
determination  of  the  reactor  volume),  and  the  volumetric  flow  rate  of  the  reactants.  Space 
velocity,  as  previously  presented  in  Equation  2.4  can  be  defined  in  several  ways.  This 
includes  gaseous  space  velocity  or  liquid  space  velocity.  Space  velocity  can  also  be 
determined  with  regards  to  a  single  specie  or  all  reactants  traveling  through  a  processor. 
Additionally,  considering  an  acoustic  field  within  a  processor,  an  effective  space  velocity 
can  be  defined  basing  an  effective  reformer  volume  on  a  typical  particle  path  length.  In 
order  to  avoid  confusion,  the  space  velocity  with  regards  to  liquid  methanol  is  used. 
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Liquid  Hourly  Space  Velocity  of  Methanol  (LHSV-M),  as  defined  in  Equation  3.4  will  be 
used  throughout  the  presentation  of  the  results. 


i.hsv-  M 


f       3/ 

/hr  ''4U'^  methanol  input 


(3.4) 


Removal  of  Data  Points 


As  multiple  sets  of  data  were  used  in  the  present  study  with  variability  in  both 
input  and  output  parameters,  output  data  were  analyzed  to  identify  outlying  points  using  a 
first  application  of  Chavenets  criterion  (Holman,1994). 

Procedures 

Specific  procedures  and  checklists  for  start-up,  data  collection,  and  shutdown 
were  made  in  order  to  collect  all  pertinent  data  as  well  as  to  avoid  hazardous  situations 
and  damage  to  equipment.  These  procedures  were  made  by  both  consulting  past  studies 
(Sterchi,  2001)  foresight  in  the  operating  process,  and  experience.  The  checklists  used, 
which  outline  the  general  procedures  for  start-up,  data  collection,  and  shutdown  can  be 
found  in  Appendix  B. 

Modeling  the  Acoustically  Enhanced  Reformation  Process 

As  discussed  in  the  theoretical  approach,  there  are  several  analytical  models  that 
may  be  applied  to  the  various  steps  of  the  reformation  process  with  acoustics.  These  are 
increase  of  effective  space  velocity  (Equation  2.1 1),  heat  transfer  (using  Equations  2.12 
and  2.13  for  Nusselt  number),  mass  transfer  (using  Equation  2.12  and  2.13  for  Sherwood 
number)  and  chemical  kinetics  (using  Equation  2.1 1  for  the  effective  area).  While  each 
of  these  models  may  represent  a  portion  of  the  overall  enhancement  due  to  the  acoustic 
waves,  an  actual  reactor  does  not  fit  the  ideal  models  developed  and  other  reactor 
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parameters  must  be  considered.  This  is  especially  true  in  the  packed  bed  reactor  utilized 
for  this  investigation  because  the  models  developed  do  not  account  for  obstruction  in  the 
acoustic  path  or  the  volume  that  the  catalyst  pellets  occupy.  While  a  comprehensive 
model  for  the  actual  conditions  may  not  be  currently  possible,  empirical  modifications  of 
the  ideal  model  for  the  experimental  facility  can  be  developed.  One  should  note  that  the 
facility  is  not  an  optimal  design  for  employment  of  acoustics  or  reforming  and  that  the 
resulting  model  is  specific  to  the  facility  used.  Although  this  is  a  significant  limitation, 
one  gains  understanding  of  how  acoustics  could  be  applied  in  different  reactors  and  just 
what  the  demonstrated  benefits  are  in  a  particular  application. 
Factorial  Experiment  Design 

A  fraction  of  the  data  collected  was  used  as  a  standard  2k  factorial  experiment 
design.  This  type  of  experiment  design  is  common  for  statistically  identifying  effects  of 
certain  factors  on  a  process  and  yields  a  model  of  the  process  with  respect  to  the  variables 
investigated  (Lawson  and  Erjavec,  2000).  A  brief  overview  of  factorial  experiment 
design  and  analysis  is  given  in  Appendix  C. 

This  factorial  experiment  design  is  used  to  develop  an  empirical  model. 
Although  the  model  developed  yields  linear  relationships  for  the  variables  investigated, 
residuals  or  curvature  of  the  output  with  relation  to  the  inputs  was  also  investigated  with 
further  data.  The  output  (dependent  variable)  for  the  factorial  experiment  design  was 
conversion  of  methanol.  After  consulting  previous  experimentation  with  the  facility, 
analysis  of  preliminary  data,  and  alignment  with  the  objectives  of  the  study,  the  input  (or 
independent)  variables  were  chosen  to  be  fuel  (premix)  How  rate,  acoustic  pressure,  and 
catalyst  bed  length.  Other  possible  variables  were  held  constant  for  the  experiments 
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including  control  setpoint  temperature  (250 °C).  frequency  (650  Hz),  and  water-methanol 
stoichiometry  (1.5:1  on  a  molar  basis).  High  and  low  levels  were  determined  for  each  of 
the  varying  parameters.  The  levels  for  fuel  flow  rate  were  a  high  of  9.5  g/min.  and  a  low 
of  4.8  g/min.  The  levels  for  acoustic  pressure  were  a  high  of  5000  Pa.  and  a  low  of  0  Pa. 
The  levels  for  catalyst  bed  length  were  a  high  of  0.576  m,  and  a  low  of  0.279  m.  This 
allowed  for  a  23  factorial  experiment  design  yielding  8  combinations  of  the  input 
variables  at  high  and  low  settings.  The  input  matrix  and  run  orders  can  be  found  for  the 
factorial  experiments  in  Table  3.5. 

Table  3.5:  Factorial  input  matrix  and  run  order. 


Flow 

Acoustic 
Pressure 

Catalyst  Bed 
Length 

Run  Order 

+ 

+ 

+ 

3,6,10 

- 

+ 

+ 

1.4.11 

+ 

- 

+ 

2,7,9 

- 

- 

+ 

5,8,12 

+ 

+ 

- 

14,18,19 

- 

+ 

- 

13,17,22 

+ 

- 

- 

15,23,24 

- 

- 

- 

16.20,21 

Although  other  parameters  were  held  "constant",  variability  in  these  settings  of 
high  and  low  levels  was  apparent.  These  fluctuations  yielded  variability  in  the  measured 
output  of  conversion.  Two  replications  of  the  study  were  performed  to  allow  statistical 
analysis  and  identification  of  outlying  points.    As  shown  by  the  run  order,  the  variables 
of  fuel  (premix)  How  rate  and  acoustic  pressure  were  randomized  while  catalyst  bed 
length  was  taken  sequentially  at  a  high  then  low  level.  The  sequential  order  of  catalyst 
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bed  length  was  required  because  of  the  time-intensive  reduction  process  between 
changing  bed  lengths.  Inadvertently  the  model  includes  variability  from  both  the 
reduction  process  and  catalyst  arrangement  as  a  bed.  This  came  as  a  result  of  rupturing 
the  seals  of  the  facility  halfway  through  the  full  (high  level)  catalyst  bed  experiments. 
The  consequences  of  this  disturbance  are  that  the  variation  in  fuel  conversion  may  be 
higher  in  the  empirical  model  than  would  be  observed  normally. 
Further  Data  and  Verification  of  the  Model 

Further  data  were  required  in  order  to  verify  the  empirical  model  constructed  and 
to  demonstrate  different  effects  of  the  acoustic  waves  on  the  reformation  process. 
Curvature  for  the  input  variables  of  the  model  was  found  by  investigating  intermediate 
sound  pressure  levels,  flow  rates  and  various  bed  lengths  including  those  bed  lengths 
lower  than  the  low  level  in  the  factorial  experiment.  This  bed  length  data  enabled  one  to 
not  only  find  the  curvature  and  limits  of  the  proposed  model  but  also  to  estimate  how  the 
local  conversion  of  methanol  progressed  though  the  full  catalyst  bed. 

Other  data  not  directly  accounted  for  in  the  empirical  model  were  also  collected. 
As  frequency  was  suspected  to  change  the  boundary  layer  thickness  for  laminar  flow,  and 
possibly  induce  the  flow  into  the  turbulent  regime,  conversion  at  two  frequencies  was 
investigated.  Because  the  temperature  within  the  processor  was  expected  to  significantly 
affect  the  conversion  of  fuel,  temperature  data  were  collected  for  acoustic  field  and  flow 
conditions.  While  holding  the  temperature  control  setpoint  at  250°C,  actual  temperature 
data  at  certain  locations  throughout  the  catalyst  bed  were  collected.  These  temperature 
data  were  collected  by  utilizing  thermocouples  extending  into  the  bed  at  specified  lengths 
through  the  unused  sampling  ports  shown  in  Figure  3.6.  Other  parameters  of  interest. 
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including  conversion  change  with  time  (indicating  catalyst  degradation),  hydrogen  power 
output,  and  acoustic  driver  input  electrical  power  were  also  collected.  Specific 
configurations  of  all  controllable  parameters  for  these  further  data  is  explained  in  the 
discussion  of  those  data  in  Chapter  6.  The  run  orders  for  the  additional  data  collected  is 
found  in  Appendix  D. 


CHAPTER  4 
PROPERTIES  OF  THE  ACOUSTIC  FIELD 

This  chapter  is  a  presentation  of  the  results  found  by  investigating  the  acoustic 
field  within  the  catalyst  bed  and  the  analysis  of  those  results.  As  the  acoustic  velocity 
must  be  higher  than  the  average  velocity  in  order  to  enhance  the  reforming  process, 
standing  acoustic  waves  built  at  resonance  with  their  associated  pressures  and  velocities 
can  be  utilized.  The  analysis  of  the  acoustic  field  provides  the  basis  for  understanding 
what  resonant  acoustic  modes  can  be  employed  in  enhancing  the  reformation  process. 
Transfer  functions  were  gathered  for  frequency  analysis.  These  data  were  compared  to 
theory  and  enabled  an  acoustic  modal  analysis  to  be  performed. 

The  Use  of  Resonance  within  the  Catalyst  Bed 

As  explained  in  the  theoretical  approach,  enhancement  of  the  reforming  process  is 
expected  from  increased  mixing  within  the  catalyst  bed  when  the  acoustic  velocity  is 
sufficiently  high.  High  amplitude  sound  waves  are  made  possible  by  building  up  standing 
or  resonant  waves  within  the  reactor.  Resonant  frequencies  and  sound  pressure  levels  in  a 
wave  guide  such  as  the  catalyst  bed  is  a  function  of  acoustic  boundary  conditions,  the 
speed  of  sound  in  the  wave  guide  and  the  damping  characteristics  of  any  obstruction 
within  the  waveguide.  A  simple  example  of  a  standing  wave  pattern  of  an  acoustic  field 
within  a  wave  guide  and  the  analogous  resonant  motion  of  a  string  is  shown  in  Figure  4. 1 . 
The  boundary  conditions  shown  in  these  simple  examples  are  the  same  as  those  expected 
within  the  catalyst  bed  studied.  The  acoustically  closed  condition  was  expected  at  the 
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Figure  4.1 :  Standing  wave  pattern  of  displacement  and  velocity  in  a  waveguide  and 
an  analogous  pattern  of  displacement  and  velocity  in  a  driven  string. 

entrance  to  the  catalyst  bed.  This  entrance  is  where  a  0.635  cm  (0.25  in.)  O.D.  tube  opens 
to  a  2.68  cm  O.  D.  (0.75  in.  schedule  40)  pipe  previously  shown  in  Figure  3.5.  Acoustic 
waves  propagate  upward  from  the  catalyst  bed  to  this  junction.  By  investigating  the 
theoretical  acoustic  power  transmitted  from  a  large  pipe  opening  to  a  small  pipe  opening 
as  shown  in  Figure  4.2,  the  expectation  of  a  closed  acoustic  condition  at  the  catalyst  bed 
entrance  can  be  justified.  The  transmitted  acoustic  power  is  shown  in  Equation  4.1, 
where  Tp  is  the  transmitted  acoustic  power,  S,  is  the  large  pipe  cross-sectional  area  and  S2 
in  the  small  pipe  cross-sectional  area  (Kinsler  et  al.,  1982).  For  the  pipe  junction 
described  above,  the  acoustic  power  transmitted  through  the  entrance  is  calculated  to  be 
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Figure  4.2:  Transmission  and  reflection  of  an  acoustic  plane  wave  at  a  pipe 

junction. 
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less  than  18%  of  the  incident  power.  Thus  acoustically  the  entrance  location  acts  as  a 
closed  end.  The  other  end  of  the  waveguide  is  where  the  acoustic  driver  is  attached. 
Although  physically  this  is  not  a  straight  through  path,  acoustically  it  appears  to  be  so.  It 
is  expected  that  because  of  the  displacement  associated  with  the  voice  coil  of  the  driver 
that  this  location  would  act  as  an  acoustically  open  condition  with  associated  large 
acoustic  displacements  and  velocities.  There  are  several  resonant  modes  and  frequencies 
associated  with  these  boundary  conditions.  The  predicted  resonant  frequencies  for  these 
closed-open  conditions  are  shown  in  Equation  4.2  where  n  is  the  mode,  c  is  sonic  velocity 
and  L  is  the  length  of  the  tube.  The  pressure  and  velocity  equations  for  any  location  and 
time  for  the  closed-open  acoustic  boundary  conditions  in  an  idealized  wave  guide  are 
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(2n-  l)c 
f,=       4L      ,     n=  1.2.3,...  (4.2) 

given  in  Equations  4.3  and  4.4  where  x  is  the  position  from  the  closed  end  (entrance)  and 
t  is  time.    With  the  theory  established  one  can  measure  the  actual  resonant  properties  of 
the  catalyst  bed  via  transfer  functions  and  perform  a  modal  analysis  to  verify'  the  predicted 
boundary  conditions. 

(2«-  \)7Dc]         I  \ 

P'„  =  Pn  cos[       2L      J  sm(2rf„t)  (4.3) 


,  .  f(2ii-l)«| 

w.  =  (7„  sin 


n  n 


COSl 


(2</)  (4.4) 


V       2L 
Transfer  Functions 
In  order  to  observe  what  resonant  frequencies  were  available  in  the  catalyst  bed 
several  transfer  functions  were  measured  in  the  reactor.  These  transfer  functions  were 
measured  from  the  signal  generator  or  driver  input  signal  to  the  pressure  transducer 
output  signal.  Excitation  of  a  range  of  frequencies  was  accomplished  by  a  swept  sine 
method.  Frequencies  were  input  from  10  Hz  up  to  15  kHz.  The  spectrums  of  the  input 
and  output  were  measured  and  compared.  This  resulted  in  a  transfer  function  magnitude 
plot.  The  transfer  function  phase  plot  and  coherence  plot  were  also  measured.  The 
coherence  was  checked  to  verify  that  the  output  signal  was  due  to  the  input  signal  and  not 
due  to  noise.  The  number  of  samples  acquired  for  the  transfer  function  was  8192.  These 
samples  were  taken  at  a  sampling  rate  of  35  ksamples/sec  resulting  in  a  frequency 
resolution  of  4.27  Hz.  Thirty  single  transfer  functions  were  averaged  to  ensure  the 
excitation  and  capture  of  all  desired  frequencies.  In  this  manner,  transfer  functions  were 
measured  for  a  variety  of  temperatures  and  reactor  conditions.  These  conditions  were 
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ambient  (-33  °C)  and  operating  temperature  (250°C  Control  Setpoint),  empty  and  full  of 
catalyst,  as  well  both  stagnant  (no  premix  flow)  and  while  generating  hydrogen  (4.9  g/min 
premix  flow)(see  Figures  4.3-  4.6).    The  transfer  functions  for  these  conditions  presented 
in  Figures  4.3  -  4.6  are  a  composite  of  two  transfer  functions  for  each  figure.  One  transfer 
functions  is  for  the  frequencies  from  40  to  2100  Hz  and  the  second  transfer  function  is  for 
2100  Hz  to  15  kHz.  There  are  no  significant  frequency  peaks  beyond  7000  Hz  in  all 
cases.  The  most  basic  of  the  configurations  tested  was  an  air  filled  ambient  temperature 
reactor  with  no  catalyst.  The  transfer  function  for  this  condition  is  found  in  Figure  4.3. 
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Figure  4.3:  Transfer  function  magnitude  for  the  catalyst  bed  housing  containing  only 

33°Cair. 

This  transfer  function  shows  the  frequencies  at  which  resonance  is  theoretically  possible. 
Transfer  functions  were  also  found  for  a  33  °C  reactor  containing  catalyst  and  air  as 
shown  in  Figure  4.4.  The  amplitude  peaks  are  lower  than  without  the  catalyst  because  of 
the  acoustic  attenuation  of  the  obstructing  catalyst  particles,  yet  the  frequency  peaks 
corresponding  to  the  highest  amplitude  modes  were  still  well  established.  After  reducing 
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Figure  4.4:  Transfer  function  magnitude  for  a  catalyst  bed  housing  containing  catalyst 

and  air  at  33  °C 

the  catalyst  and  holding  carbon  dioxide  in  the  reactor  at  250 °C,  the  transfer  function,  as 
shown  in  Figure  4.5  revealed  a  small  shift  in  frequencies  due  to  a  slight  change  of  the 
sonic  velocity  of  the  gas  within  the  catalyst  bed.  The  sonic  velocity  in  33  °C  air  is  351  m/s 
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igurc  4.5:  Transfer  function  magnitude  for  a  reactor  containing  reduced  catalyst  and 
C02  at  operating  temperature  (250  °C). 
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while  the  sonic  velocity  in  250°C  carbon  dioxide  is  357  m/s.  With  this  slight  difference 
in  sonic  velocity  only  a  small  change  in  frequency  positioning  is  expected.  These 
frequency  peaks,  under  the  conditions  described  were  also  found  positioned  at  similar 
relative  locations  as  the  air  filled  reactor  with  catalyst,  indicating  no  shift  in  primary 
oscillation  modes.  With  a  flow  of  reactants  and  products  in  the  catalyst  bed  (premix  flow 
of  4.8  g/min)  under  operating  conditions  (250°C  control  setpoint)  Transfer  functions 
were  acquired  as  displayed  in  Figure  4.6.  As  seen  in  this  figure,  the  transfer  function 
shifted  slightly  in  magnitude  and  frequency  positioning  when  water/methanol  premix  was 
being  reacted  because  of  both  the  slight  differences  in  overall  reactor  temperature, 
pressure,  and  gas  composition  (and  thus  sonic  velocity)  due  to  the  reaction.  The  highest 
amplitude  peak  previously  near  550  Hz  was  shifted  to  650  Hz  and  the  peak  previously  at 
1700  Hz  shifted  to  1900  Hz,  while  two  low  amplitude  peaks  became  evident  near  1400 
Hz  and  2300  Hz  .  An  example  of  the  time  domain  pressure  signal  within  the  catalyst  bed 
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Figure  4.6:  Reactor  transfer  function  magnitude  while  generating  hydrogen  at  250 °C. 
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while  being  driven  at  650  Hz  and  while  reaction  is  taking  place  was  presented  previously 
in  Figure  3.14.  It  was  determined  that  while  generating  hydrogen  rich  products,  an 
acoustic  sound  pressure  level  of  165  dB  (re20  uPa)  could  be  maintained  in  the  catalyst 
bed  if  driven  at  650  Hz  and  approximately  9  watts  of  electrical  power  to  the  driver. 

One  can  compare  these  experimental  results  of  resonant  frequencies  found  by  the 
transfer  functions  with  theoretical  resonant  frequencies  for  the  assumed  boundary 
conditions.  Allowable  resonant  frequencies  (f)  for  an  closed-open  waveguide  have  been 
given  in  Equation  4.1.  Noticeably,  the  frequencies  as  predicted  by  Equation  4.1  are  not 
all  present  in  the  measured  transfer  functions.  This  deviation  from  the  analytical  solution 
is  due  to  insufficient  driving  power  at  certain  frequencies  and/or  to  the  actual  boundary 
conditions  present  in  the  reactor.  The  midrange  driver  has  a  frequency  response  as  shown 
in  Figure  4.7  (Klipsch,  2001).  Because  the  lower  frequencies  (below  200  Hz)  are  not 
well  excited  by  the  midrange  horn  driver,  it  was  not  expected  that  the  lower  frequencies 
would  be  excited  in  the  reactor. 
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Figure  4.7:  Frequency  response  of  the  acoustic  driver  as  measured  in  a  plane 
wave  tube  by  the  manufacturer. 
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The  actual  resonant  frequencies  from  the  transfer  functions  can  be  translated  to 
their  associated  wavelengths  (X)  through  Equation  4.2  With  the  wavelengths  and  at  least 
one  boundary  condition  known,  modal  shapes  can  be  plotted  and  the  location  of  actual 
boundary  conditions  can  be  found. 


J  n 

Modal  Analysis 


(4.2) 


Wavelengths  in  the  33°C  air-filled  catalyst  bed  housing  without  catalyst  were 
calculated  from  the  measured  transfer  function  previously  shown  in  Figure  4.3.  With  the 
acoustically  closed  condition  assumed  as  the  catalyst  bed  housing  entrance,  the  velocity 
and  displacement  waveforms  for  six  frequencies  (amplitude  has  been  normalized)  were 
plotted.  These  plotted  waveforms  are  shown  in  Figure  4.8.  This  figure  plots  normalized 
velocity  magnitude  versus  position  from  the  entrance,  and  has  the  temperature  and 
pressure  ports  along  with  suspected  physical  boundary  conditions  marked  along  the  x- 
axis. 
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Figure  4.8:  Normalized  velocity  magnitude  for  six  waves  found  in  a  33°C  reactor  without 

catalyst. 
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From  Figure  4.8,  it  is  evident  that  an  acoustically  open  condition  exists  near 
800-810  mm  inside  the  reactor.  This  open  condition  directly  corresponds  to  the  driver 
mounting  plate  where  the  driver  is  attached  to  the  pipes  leading  to  the  catalyst  bed.  This 
location  is  termed  the  adapter  throat.  This  acoustically  open  condition  at  this  point 
indicates  that  the  mass  of  fluid  inside  the  acoustic  driver  oscillates  in  unison  with  the 
voice  coils  and  diaphragm.  Considering  the  large  diameter  (7.62  cm)  of  the  diaphragm 
and  the  relatively  short  distance  from  the  diaphragm  to  the  adapter  throat  (8.9  cm),  this 
unison  in  oscillation  is  not  surprising.  In  addition  to  this  boundary  condition  at  the 
adapter  throat,  there  are  other  contributing  boundary  conditions  for  some  of  the  frequency 
peaks  found  in  the  transfer  functions.  Most  notably  at  550  Hz,  1 700  Hz,  and  2100  Hz  for 
the  empty  air  filled  reactor  there  is  contribution  from  an  open  condition  at  the  pressure 
transducer  (PT)  tube  and  housing  (see  Figure  3.6).  As  expected,  the  frequency  peaks  with 
contributing  boundary  conditions  are  higher  than  those  without.  This  amplification  of 
certain  frequencies  with  contributing  boundary  conditions  explains  why  the  real  reactor 
deviates  from  the  idealized  case  presented  by  Equation  4.2.    As  can  be  seen  from  Figure 
4.6,  the  contributing  conditions  also  widen  the  frequency  peaks  which  allows  for  better 
control  of  the  sound  pressure  level.  As  noted  in  the  explanation  of  the  experimental 
facility  there  are  fluctuations  in  overall  temperature  and  gas  composition  with  time  that 
can  change  the  resonant  characteristics  of  the  catalyst  bed.  The  wide  bands  of  resonance 
found  for  frequencies  with  contributing  boundary  conditions  allow  for  minimization  of 
the  changes  in  sound  pressure  when  acoustically  driving  at  these  frequencies. 

The  modal  analysis  was  verified  by  moving  the  pressure  transducer  down  the 
reactor  length  through  the  six  access  ports  in  the  air-filled  reactor  without  catalyst.    The 
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transfer  functions,  measured  through  the  ports  are  shown  in  Appendix  E  for  the  lower 
frequencies.  As  expected  from  the  standing  wave  pressure  amplitude  trace,  which  can  be 
thought  of  as  90°  leading  the  velocity  trace  in  Figure  4.8,  the  transducer  picked  up  higher 
magnitudes  pressures  at  the  corresponding  pressure  anti-nodes.  These  shifts  in  the 
transfer  function  magnitude  and  also  in  the  corresponding  phase  plots  were  noted  and 
found  to  follow  the  expected  modes  as  predicted  by  the  waveforms  in  Figure  4.8.    A 
summary  of  the  frequencies,  wavelengths  and  boundary  conditions  for  the  air-filled  33  °C 
reactor  without  catalyst  is  shown  in  Table  4. 1  below. 
Table  4.1 :  Waveforms  and  boundary  conditions  for  a  33  °C  reactor  containing  only  air. 
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0.11 
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0.10 
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0.09 
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0.09 

Open  at  driver  throat 

Open  at  PT  Tube 

5383 

0.06 

Open  at  adapter  throat 

Open  at  PT  Tube 

6669 

0.05 

Closed-closed  transverse  wave 

Open  at  PT  Tube 

velocity  amplitude  mode  shape  is  plotted  in  Figure  4.8 


CHAPTER  5 
EMPIRICAL  MODEL  DEVELOPMENT  AND  ANALYSIS 

This  chapter  presents  and  discusses  the  results  of  the  23  factorial  experiment 
design  used  to  develop  empirical  correlation  factors  to  model  the  reformation  process 
under  the  influence  of  a  resonant  acoustic  field.  The  reactor  was  investigated  for  the 
input  variables  of  premix  flow  rate,  sound  pressure  amplitude,  and  catalyst  bed  length  at 
high  and  low  levels  as  described  in  the  experimental  approach.  This  investigation  allowed 
for  correlation  factors  of  an  empirical  model  to  be  found.  Model  correlation  factors  are 
presented,  along  with  an  uncertainty  analysis  of  those  correlation  factors.  Midpoint 
curvature  for  each  variable  and  the  effect  of  frequency  on  the  conversion  level  of  the 
steam-reforming  process  were  also  examined. 

Experimental  Results  and  Model  Correlation  Factors 

By  using  a  standard  factorial  design  of  experiments  one  can  develop  an  empirical 
model  of  a  certain  design  space  with  relatively  few  experimental  trials.  An  overview  of 
the  general  setup  of  such  a  process,  formulation  of  the  correlation  factors,  and  statistical 
analysis  of  hypothetical  results  are  presented  in  Appendix  C.  Actual  values  for  the  high 
(+)  and  low  (-)  levels  investigated  and  the  randomized  run  order  were  previously 
presented  in  Chapter  3.  The  results  of  the  experiments  for  the  configurations  specified 
yielded  the  output  data  shown  in  Table  5.1.  In  these  experiments  the  output  or 
independent  variable  is  conversion  of  methanol  in  percent.  The  temperature  control 
setpoint  was  held  constant  as  were  all  other  controllable  parameters  not  specified.  As 
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shown  in  Table  5.1  the  standard  deviations  are  relatively  low  indicating  good 

repeatability  for  each  configuration.  The  data  also  show  that  with  only  a  few  exceptions 

the  conversion  decreases  over  the  trials  run.  Because  the  trial  run  is  correlated  to  the 

relative  sample  time,  these  data  suggest  that  catalyst  degradation  was  causing  the  bulk  of 

variability. 

Table  5.1 :  The  input  matrix  and  methanol  conversion  percentage  data  obtained  from  the 

factorial  experiments. 


Run 

Flow 
Rate 

(XI) 

Sound 
Pressure 

(X2) 

Catalyst 

Bed 
Length 

(X3) 

Trial 
1 

Trial 

2 

Trial 

3 

Conversion 

%  Average 

(Y) 

Standard 
Deviation 

1 

+ 

+ 

+ 

96.33 

96.28 

94.94 

95.85 

0.79 

2 

- 

+ 

+ 

99.74 

99.61 

99.44 

99.60 

0.15 

3 

+ 

- 

+ 

93.06 

92.22 

91.84 

92.37 

0.62 

4 

- 

- 

+ 

99.13 

98.81 

98.90 

98.95 

0.17 

5 

+ 

+ 

- 

70.58 

70.19 

70.34 

70.37 

0.20 

6 

- 

+ 

- 

92.57 

91.46 

91.33 

91.79 

0.68 

7 

+ 

- 

- 

61.74 

60.82 

59.77 

60.78 

0.99 

8 

- 

- 

- 

87.61 

86.45 

86.37 

86.81 

0.69 

Equation  5.1  is  used  to  find  the  standard  error  (SE)  from  the  pooled  standard 
deviation  (op)  and  the  total  number  of  factorials  (nf)  for  these  data.  The  total  average 

conversion  (  Y  ),  pooled  standard  deviation,  standard  error,  and  total  degrees  of  freedom 

for  the  data  are  given  in  Table  5.2.    These  data  suggest  that  mean  values  within  the 
conditions  studied  are  expected  to  deviate  less  than  or  equal  to  the  standard  error  (0.251 
%)  and  any  individual  point  should  not  deviate  by  more  than  twice  the  pooled  standard 
deviation  (1.22  %)  from  its  mean  value  using  a  95%  confidence  interval.  The  correlation 


sE  = 


2a, 
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(5.1) 


Table  5.2:  Total  average  conversion,  pooled  standard  deviation,  standard  error,  and  total 
degrees  of  freedom  from  the  factorial  experiments. 

Total  Average 
Conversion  % 

Pooled  Standard  Deviation 

Standard  Error 

Total  Degrees  of 
Freedom 

87.06 

0.61 

0.251 

16 

factors  or  effects  (E)  and  interactions  (I)  were  calculated  for  each  variable  (X)  and  set  of 
variables  by  taking  the  average  of  the  average  conversion  outputs  (Y)  when  the  variable 
or  set  of  variables  were  high  and  subtracting  the  same  average  when  the  variable  or  set  of 
variables  were  low.  This  statement  can  be  simplified  to  yield  Equation  5.2  for  the  effects. 


rp  Run=\ 


(5.2) 


Interactions  are  found  in  a  similar  manner  where  the  combination  of  the  level  X,  run 
multiplied  by  X,  run  is  used  instead  of  X,  run  in  Equation  5.2  find  the  Interaction  (L).  This 
calculation,  using  the  data  presented  in  Table  5.1  resulted  in  the  correlation  factors  for  the 
model  as  the  effects  and  interactions  displayed  in  Table  5.3 

Table  5.3:  The  effects  and  interactions  found  from  the  factorial  experimentation. 


Effect  of  Flow  Rate  (E,) 

-14.44 

Effect  of  Sound  Pressure  (E2) 

4.67 

Effect  of  Catalyst  Bed  Length  (E3) 

19.26 

Interaction  of  Flow  Rate  and  Sound  Pressure  (I,-,) 

1.86 

Interaction  of  Flow  rate  and  Catalyst  Bed  Length  (I13) 

9.28 

Interaction  of  Sound  Pressure  and  Catalyst  Bed  Length  (I23) 

-2.61 

Interaction  of  Flow  Rate,  Sound  Pressure  and  Catalyst  Bed  Length  (Il23) 

-0.45 
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These  correlation  factors  show  the  magnitude  and  direction  each  variable  has  on 

the  conversion  of  methanol  in  the  reformer.  These  model  correlation  factors  yield  the 

equation  for  methanol  conversion  (C)  as  shown  in  Equation  5.3.  However,  not  all  effects 

F  E  E  I  I  I  1 

C=  Y  +  ^-X,  +  ^fX,  +  ^rX,  +  -^X,X.  +  -^X,X,  +  -^X2X3  +  -J2-  X,X2X,(53) 

2222222 

or  interactions  presented  above  are  necessarily  statistically  significant  or  need  to  be 
included  in  the  working  model  equation.  In  order  to  find  the  statistical  significance  of  the 
correlation  factors,  one  must  perform  an  uncertainty  analysis  on  the  effects  and 
interactions. 

Uncertainty  Analysis 
In  order  to  formulate  the  proper  model  correlation  equation  one  must  test  each 
effect  and  interaction  for  statistical  significance.  This  uncertainty  analysis  can  also  be 
used  to  test  the  significance  level  of  the  individual  correlation  factors.  This  is  particularly 
important  for  statistically  proving  that  the  acoustic  wave  has  a  significant  effect  on  the 
conversion  in  this  reactor.  The  signal  to  noise  t-ratio  (t*)  was  found  using  Equation  5.4., 

E  or  I 
t*  -  —  (5.4) 

As  with  a  standard  t-test,  statistical  significance  for  each  effect  and  interaction  was  found 
at  the  desired  confidence  level  (both  95%  and  99.9%  are  presented)  for  the  total  degrees 
of  freedom,  by  comparing  the  absolute  value  of  the  signal-to-noise  t-ratio  to  the  Student  t- 
value.  The  Student-t  value  yields  the  maximum  noise  ratio  at  the  specified  degrees  of 
freedom  and  confidence  interval.  The  effects  and  interactions  with  absolute  values  of  the 
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signal  to  noise  t-ratios  larger  than  the  Student  t-value  are  considered  to  be  significant  at 
the  confidence  interval  stated.  The  results  of  this  analysis  are  shown  in  Table  5.4. 
Table  5.4:  Signal-to-noise  t-ratios  and  statistical  significance. 


t* 

Significant  at 
95%  Confidence? 

Significant  at  99.9% 
Confidence? 

Student  t-value  at  95%  Confidence 

2.120 

Student  t-value  at  99.9%  Confidence 

4.015 

Effect  of  Flow  Rate  (E,) 

-57.60 

Yes 

Yes 

Effect  of  Sound  Pressure  (E2) 

18.64 

Yes 

Yes 

Effect  of  Catalyst  Bed  Length  (E3) 

76.80 

Yes 

Yes 

Interaction  of  Flow  Rate  and  Sound 
Pressure  (I,2) 

7.42 

Yes 

Yes 

Interaction  of  Flow  Rate  and  Catalyst 
Bed  Length  (I13) 

37.02 

Yes 

Yes 

Interaction  of  Sound  Pressure  and 
Catalyst  Bed  Length  (I23) 

-10.41 

Yes 

Yes 

Interaction  of  Flow  Rate,  Sound 
Pressure  and  Catalyst  Bed  Length  (I123) 

-1.78 

No 

No 

These  data  also  show  that  all  of  the  variables  investigated  including  sound 
pressure  has  a  statistically  significant  effect  on  the  conversion  for  this  reactor  at  a 
confidence  level  of  99.9%.  The  fact  that  the  interaction  of  flow  rate,  sound  pressure  and 
catalyst  bed  length  can  be  neglected  allows  for  a  slight  simplification  of  the  conversion 
equation.  The  statistically  valid  conversion  model  for  the  inputs  of  flow  rate,  sound 
pressure  amplitude  and  catalyst  bed  length  is  found  in  Equation  5.5.  This  equation  is 
used  by  finding  where  the  actual  flow  rate,  sound  pressure  and  bed  length  variables  lit 
into  the  high  (+1 )  and  low  (-1 )  settings  for  each  parameter.    These  values  can  be  found  by 


C(%)  =  87.06  -  7.22X,  +  2.34X:  +  9.63X, 
+  0.93X.X,  +  4.64X,X,  -  1.31X2X3 


(5.5) 
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utilizing  Equations  5.6  -  5.8  for  specific  premix  flow  rates  (g/min),  sound  pressure 
amplitudes  (Pa),  and  catalyst  bed  lengths  (m),  respectively.  These  values  are  then  utilized 
in  Equation  5.5  as  X„  X2,  and  X3  where  C  is  the  conversion  of  methanol  in  percent. 

_  Flow(g  /  min)-  7.13(g  /  min) 
1  "  2.35(g/min) 

Pressure(Pa)-2680(Pa) 
X:  "  "  2680(Pa)  (5*7) 

Bed  Length  (m)-  0.424  (m) 

x3 " fsioio <5-8> 

While  sound  pressure  (X2)  does  not  influence  the  conversion  as  much  as  the  flow  rates 
(X,)  or  catalyst  bed  lengths  (X3)  studied,  the  equation  does  show  how  much  an  effect 
sound  pressure  can  have  in  relation  to  the  other  variables.  In  order  to  demonstrate  the 
some  of  the  limitations  of  the  model,  midpoint  curvature  data  were  obtained  and  the 
effect  of  frequency  was  observed. 

Midpoint  Curvature 

The  model  presented  in  Equation  5.5  was  further  examined  for  curvature  near  the 
center  points  of  the  variables  investigated. 
Flow  Rate 

Flow  rate  curvature  in  the  model  was  investigated  at  two  sound  pressure  levels 
keeping  catalyst  bed  length  constant  as  shown  in  Figure  5.1 .  For  these  investigations, 
catalyst  bed  length  was  set  at  0.568  m  and  the  high  level  of  sound  pressure  was 
approximately  5000  Pa  or  165  dB.  Four  intermediate  points  were  taken  between  the  high 
flow  setting  (+1)  of  9.5  g/min  and  the  low  flow  setting  (-1)  of  4.8  g/min.  Comparing  the 
conversion  at  various  flow  rates  with  the  conversion  predicted  by  the  model,  for  cases 
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Figure  5.1 :  Predicted  and  actual  conversions  for  various  flow  rates  with  and  without 

acoustics. 

with  and  without  the  1 65dB  acoustic  field  show  some  visual  curvature  near  the 
midpoints.  Although  one  can  see  this  slight  trend  upward  in  conversion  at  the  midpoints, 
no  statistically  significant  curvature  is  found.  This  is  because  differences  between  the 
predicted  and  actual  conversion  data  (also  known  as  the  residuals)  are  in  the  acceptable 
noise  range  (±  1 .22  %)  for  the  individual  data  points  collected. 
Sound  Pressure 

Curvature  for  the  sound  pressure  was  investigated  at  a  high  level  of  premix  How 
rate  (9.5  g/min)  and  a  high  level  of  catalyst  bed  length  (0.576  m)  at  a  frequency  of  650 
Hz.  The  actual  conversion  at  various  sound  pressure  amplitudes,  together  with  the 
conversion  predicted  by  the  model,  are  shown  in  Figure  5.2.  Although  some  of  the  actual 
values  are  slightly  lower  than  expected,  these  values  do  not  show  significant  curvature  as 
the  residuals  are  less  than  1%  for  the  cases  investigated. 
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Figure  5.2:  Predicted  and  actual  conversion  for  maximum  pressure  amplitude  with  a 
0.568  m  catalyst  bed  length  and  9.5  g/min  premix  flow  rate. 

Catalyst  Bed  Length 

Curvature  in  the  catalyst  bed  length  variable  was  investigated  at  the  high  level  of 
premix  flow  rate  (9.48  g/min)  with  a  1 65  dB  (re20Pa)  acoustic  field  as  shown  in  Figure 
5.3.  This  figure  shows  the  predicted  conversion  of  the  model  and  actual  conversion  data 
points  gathered  at  different  catalyst  bed  lengths.  Note  that  the  conversion  predicted  by  the 
model  has  been  extrapolated  beyond  the  low  level  of  catalyst  bed  length  (0.278  m)  in  the 
figure.  Data  shown  has  significant  curvature  as  it  deviates  from  the  model  predictions 
away  from  the  high  (0.568  m)  and  low  test  levels  (0.278  m).  The  linear  model  prediction 
is  in  error  for  catalyst  bed  length  and  it  appears  that  a  logarithmic  fit,  shown  as  the 
regression  line  on  the  figure  is  more  suitable  and  realistic,  especially  considering  the 
actual  entrance  conditions.  Further  data  on  conversion  with  respect  to  catalyst  bed  length 
is  presented  in  Chapter  6. 
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Figure  5.3:  Predicted  and  actual  conversions  for  bed  length  at  a  flow  rate  of  9.5  g/min 

premix  and  a  1 65  dB  acoustic  field. 

The  Effect  of  Frequency 
As  previously  described  in  the  theoretical  approach,  for  an  ideal  setting,  frequency 
has  been  linked  to  boundary  layer  thickness  for  the  laminar  regime  (Equation  2.7)  and  the 
kinetic  Reynolds  number  (Equation  2.8).  Because  of  this,  the  acoustically  enhanced 
methanol-steam  reforming  process  was  also  investigated  to  find  any  correlation  between 
frequency  and  conversion  of  fuel.  Frequency  effects  were  measured  at  a  high  level  of 
premix  flow  rate  (9.5  g/min)  and  a  high  level  catalyst  bed  length  (0.568  m).  With  the 
acoustic  driver  power  held  constant  at  9.3  Watts,  the  average  observed  sound  pressure 
level  dropped  to  1434  Pa  when  the  reference  frequency  was  increased  to  1980  Hz.  High 
acoustic  pressures  (5000  Pa  Max.  Amplitudes)  were  not  possible  at  this  frequency  with 
the  same  driver  power  due  to  the  acoustic  resonant  properties  as  previously  presented. 
Conversion,  as  predicted  by  the  model  created  from  650  Hz  data  for  this  sound  pressure 
amplitude  was  93.2%.  This  value  matches  the  experimental  conversion  found  at  1980  Hz 
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exactly.  However,  care  should  be  taken  in  the  interpretation  of  this  data  point  because 

this  value  cannot  necessarily  be  considered  out  of  the  noise  band  of  the  data  without 

acoustics.  Notwithstanding,  this  data  suggests  that  frequency  was  not  a  significant 

variable  for  the  acoustic  fields  investigated. 

Reynolds  numbers  were  estimated  for  different  frequencies  using  the  properties  of 

ambient  pressure  steam  at  250°C  (v  =  4.401 8  x  10"5  m2/s)  and  the  diameter  of  a  catalyst 

pellet  (0.54  cm)  as  the  characteristic  length.    The  average  velocity  was  estimated  at  0.89 

m/s  from  the  flow  rate,  density  steam  at  250  °C,  and  volume  of  the  reactor.  The  Reynolds 

numbers  are  shown  in  Table  5.5.  Static  Reynolds  number  (Res)  for  the  non-oscillating 

flow  indicates  a  laminar  flow  regime  with  no  acoustics.  Kinetic  Reynolds  numbers 

(Rek)were  calculated  for  an  oscillating  field  at  the  two  frequencies  investigated  using  the 

properties  above  in  Equation  2.8.  The  kinetic  Reynolds  number  at  650  Hz  is  2705 

suggesting  that  the  flow  is  in  the  transitional  or  turbulent  regime  at  this  frequency.  An 

increase  of  frequency  to  1980  Hz  increases  the  kinetic  Reynolds  number  to  8241  which  is 

also  in  the  turbulent  regime.  A  comparison  of  the  turbulent  regime  to  turbulent  regime 

Table  5.5:  Static  and  kinematic  Reynolds  numbers  and  methanol  conversion  levels  for 

three  frequencies. 


Hz 

Reynolds  Number 

Conversion  % 

0 

110(Res) 

92.4 

650 

2705  (Rek) 

93.2  (predicted  from  model  at  1434  Pa) 

1980 

8241  (Rek) 

93.2  actual  (1434  Pa) 

with  the  two  frequencies  yields  no  difference  in  conversion  as  shown  by  the  data. 
Furthermore,  increased  frequency  yields  lower  pressure  for  the  same  driver  input  power 
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thus  the  overall  power  benefit  was  adversely  influenced  by  increasing  acoustic  frequency 
in  this  system. 

Summary 
Factorial  design  of  experiments  has  been  used  to  find  the  effects  of  the  parameters 
of  premix  flow  rate,  sound  pressure  amplitude,  and  catalyst  bed  length.  Sound  pressure 
amplitude  in  the  study  was  statistically  significant  at  a  99.9%  confidence  level.  From  the 
data  collected  a  model  was  formulated  that  predicts  methanol  conversion  for  the  above 
parameters.  The  model  had  negligible  curvature  for  fuel  flow  rate  and  maximum  pressure 
amplitude  and  significant  error  for  catalyst  bed  lengths.  Conversion  for  catalyst  bed 
length  is  better  modeled  by  a  logarithmic  relationship.  Acoustic  field  frequency  did  not 
produce  a  noticeable  change  in  conversion. 


CHAPTER  6 
THE  EFFECTS  OF  ACOUSTICS  ON  THE  REFORMATION  PROCESS 

The  acoustic  wave  affects  several  measurable  parameters  of  the  steam  reformer 
and  this  chapter  further  quantifies  these  effects.  The  output  parameters  presented  are 
axial  and  radial  temperature  profile,  conversion,  power  demand,  controllability,  catalyst 
degradation,  and  output  power  of  the  processor.  These  effects  translate  into  size 
reduction  and/or  increased  capacity  for  an  acoustically  enhanced  methanol  steam 
reformer,  and  they  point  to  potential  for  application  of  acoustic  enhancement  in  fuel  cell 
vehicles. 

The  Effect  on  Temperature  Profile 

The  steam  reformer  was  investigated  to  quantify  the  change  in  temperature  profile 
that  was  expected  due  to  increased  heat  and  mass  transfer  in  the  reactor  with  acoustic 
pulsations.  With  acoustics,  it  was  expected  that  the  radial  temperature  profile  would  have 
a  less  severe  gradient  and  the  measured  wall  temperature  would  more  closely  represent 
the  average  temperature.  Axially,  it  was  thought  that  the  endothermic  reforming  section 
would  compact  due  to  increased  heat  transfer  and  mixing.  Furthermore,  it  was  thought 
that  a  higher  temperature  section  would  move  into  the  middle  of  the  reactor  bed  with  a 
decrease  in  energy  consumed  at  that  point  signifying  the  exothermic  water  gas  shift 
reaction.  Temperature  profiles  for  various  space  velocities  and  acoustics  sound  levels  are 
given  in  Figures  6.1-  6.4.  Actual  data  points  and  the  linear  regressions  used  to  calculate 
the  profiles  for  these  figures  are  found  in  Appendix  F.  It  should  be  noted  that  the  profiles 
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given  in  figures  6.1-6.4  are  one-half  of  the  reactor  radially.  Temperatures  were  measured 
randomly  from  the  opposite  radial  position  and  the  mean  temperature  was  found  to  be 
within  the  uncertainty  (2°C)  for  the  positions  and  conditions  measured.  It  is  also 
important  to  note  that  the  temperature  profiles  represent  two-thirds  of  the  entire  axial 
length  of  the  reactor  as  there  are  heating  sections  both  above  and  below  the  profile  area 
displayed.  The  absence  of  available  sampling  ports  at  the  upper  and  lower  locations 
prevents  temperature  acquisition  at  these  positions.  The  profile  may  also  have  some  local 
inaccuracies  as  the  temperature  at  the  area  between  the  zone  exits  was  calculated  with  a 
linear  interpolation  between  the  exits.  Notwithstanding  these  limitations,  there  is  a 
noticeable  difference  between  the  measured  profiles  with  acoustics  and  those  without 
acoustics.  One  can  see  a  clear  increase  in  center  line  temperature  for  those  with  acoustics 
at  both  space  velocities.  In  the  axial  direction,  ones  sees  the  expected  increase  of 
temperature  with  acoustics  and  radially  the  profile  is  flattened,  thus  the  overall 
temperature  is  increased.  For  the  cases  with  acoustics  (Figures  6.2  and  6.4),  one  can  see 
that  the  heat  is  being  transferred  more  efficiently  throughout  the  reactor  displayed. 

For  Figure  6.1  in  the  axial  direction,  one  can  see  that  the  235  °C  isotherm  extends 
below  Zone  3  Exit.  This  indicates  that  there  is  an  endothermic  reaction  at  and  possibly 
below  this  point.  Figure  6.2  shows  the  235  °C  isotherm  ending  between  zone  exits  2  and 
3  for  this  case  with  acoustics.  As  the  endothermic  portion  of  steam  reforming  is  thought 
to  be  the  rate  limiting  step,  the  temperature  in  the  reactor  may  be  a  good  indication  of 
how  the  reaction  is  taking  place.    The  increase  in  temperature  at  the  center  line  with 
acoustics  reveals  that  there  may  be  a  chance  to  decrease  the  required  length  of  the  reactor. 
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Figure  6. 1 :  Temperature  profile  for  the  Figure  6.2:  Temperature  profile  for  the 

catalyst  bed  at  premix  flow  rate  4.9  g/min    catalyst  bed  at  premix  flow  rate  of  4.9  g/min 
without  an  acoustic  field.  with  an  acoustic  field  at  165  dB. 
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Figure  6.3:  Temperature  profile  for  the 

catalyst  bed  at  a  premix  flow  rate  of  9.5 

g/min  without  an  acoustic  field. 


Zone  1  Exit 


Zone  2  Exit 


Zone  3  Exit 


Zone  4  Exit 


Zone  5  Exit 


Center 


Wall 


Figure  6.4:  Temperature  profile  for  the 

catalyst  bed  at  premix  flow  rate  of  9.5  g/min 

with  an  acoustic  field  at  165  dB. 
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Similarly,  comparison  of  Figures  6.3  and  6.4  show  that  the  axial  temperature  is  markedly 
increased  by  the  acoustic  wave  at  higher  space  velocities.  In  order  to  verify  the  suspected 
increased  reaction  rate  as  manifested  by  the  temperature  profiles  one  must  observe  the 
conversion  of  fuel. 

The  Effect  on  Conversion 
As  the  acoustic  wave  increases  the  mixing  in  the  reactor,  a  combination  of 
increased  temperature,  and  increased  heat  and  mass  transfer  allows  for  more  fuel  to  be 
processed  by  the  reactor.  As  proven  by  the  analysis  and  presentation  in  the  preceding 
chapter,  the  conversion  of  fuel  is  significantly  enhanced  by  the  acoustic  wave  at  a  variety 
of  space  velocities.  The  conversion  of  the  fuel  and  thus  the  activity  of  the  catalyst  was 
expected  to  improve  permitting  an  increase  in  the  capacity  of  the  processor.  Conversion 
of  methanol  with  various  catalyst  bed  lengths  (and  corresponding  masses)  were  tested  to 
investigate  the  expected  increase.  Additionally  the  changing  bed  lengths  provided  an 
estimate  of  how  the  conversion  takes  place  throughout  a  full  catalyst  bed.  Table  6.1 
shows  the  results  of  the  conversion  of  methanol  with  different  bed  lengths.  One  can  see 
this  table  as  the  fuel  conversion  as  a  function  of  axial  position  or  bed  length,  although  this 
may  not  be  completely  accurate  in  interpretation.  As  these  data  include  the  6  data  points 
previously  shown  as  part  of  the  model  discussion,  this  table  gives  a  composite  map  that 
yields  conversion  throughout  the  catalyst  bed  for  different  acoustic  and  flow  rate 
parameters.  Each  value  is  the  average  of  3  experimental  trials  for  the  given  parameters. 
It  is  important  to  note  that  some  inaccuracies  might  exist  in  interpretation  of  this  table. 
Being  a  down-flow  reactor,  data  were  collected  as  the  bed  position  lowered  in  the  catalyst 
housing.  Thus  the  actual  conversion  taking  place  through  a  full  catalyst  bed  may  not  be 
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Fable  6.1 :  Extent  of  conversion 

for  catalyst  bed  length  with  and  without  acoustics. 

Premix  flow  rate  of  4.8  g/min 

Premix  Flow  Rate  of  9.5  g/min 

Catalyst  Bed 
Length 

With 
Acoustics 

Without 
Acoustics 

With 
Acoustics 

Without 
Acoustics 

1/6 

65.4  % 

55.9  % 

33.8  % 

23.6  % 

2/6 

78.6  % 

68.8  % 

46.5  % 

35.9  % 

3/6 

91.8% 

86.4  % 

70.4  % 

60.8  % 

4/6 

97.6  % 

94.8  % 

86.4  % 

78.2  % 

5/6 

99.2  % 

97.6  % 

93.3  % 

87.1  % 

6/6 

99.6  % 

99.0  % 

95.9% 

92.4  % 

represented  accurately  by  the  table.  Because  the  acoustic  wave  pattern  changes 
throughout  the  bed  length  as  discussed  in  the  modal  analysis,  the  top  l/6th  portion  of  the 
reactor  may  see  a  different  acoustic  pattern  than  the  lower  l/6th  portion  where  the  data 
was  collected.  This  same  problem  would  exist  for  the  non-oscillating  reactant  flow 
pattern  in  the  catalyst  bed.  Notwithstanding  these  difficulties.  Table  6.1  is  clear  evidence 
that  a  reduction  of  catalyst  bed  length  can  occur  by  employing  an  acoustic  field.  For  the 
case  of  4.8  g/min  of  premix  with  acoustics,  the  98%  conversion  is  realized  with  one 
quarter  less  reactor  volume  than  without  using  acoustics.  The  95%  level  at  9.48  g/min  of 
premix  shows  a  similar  reduction  in  the  required  volume  of  catalyst.  This  increasing 
difference  in  conversions  between  oscillating  and  non-oscillating  continues  at  lower 
acceptable  conversions  as  demonstrated  with  higher  flow  rates  or  smaller  reactor 
volumes. 

Figure  6.5  is  a  map  of  the  data  in  Table  6.1  showing  the  exponential  nature  of 
conversion  toward  100%  as  a  function  of  catalyst  bed  length.  One  notes  that  the 
difference  in  conversion  between  the  cases  with  and  without  acoustics  is  dramatic  for  the 


90 


♦  4.8  g/min  w/  Acoustics 
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Figure  6.5:  Conversion  versus  bed  length  with  and  without  acoustic  enhancement. 
1st  zone  with  little  change  in  this  difference  until  the  middle  zone  is  reached.  Upon 
examination  of  the  mode  shapes  previously  shown  in  Figure  4.8,  and  considering  catalyst 
within  the  reactor  from  Port  5  (also  referred  to  as  Zone  4  Exit)  to  the  end  of  the  reformer 
expansion  with  a  mode  shape  shown  as  Wave  2,  the  data  suggest  that  the  difference  in 
conversion  is  implemented  primarily  through  the  velocity  anti-nodes  of  the  acoustic  field. 
The  velocity  anti-nodes  of  the  acoustic  field  are  where  the  highest  oscillating  velocities 
and  corresponding  displacements  are  found.  These  high  velocities  are  the  primary 
mechanism  by  which  the  mixing  is  enhanced  within  the  oscillating  flow  reactor.  A  plot  of 
the  difference  in  conversion  for  cases  with  acoustic  fields  minus  the  cases  without 
acoustic  fields  for  both  flow  rates  is  shown  in  Figure  6.6  as  bed  length  increases.  The 
effect  of  the  acoustic  mode  shape  diminishes  as  the  bed  length  increases  and  conversion 
approaches  the  limiting  value  of  100%. 

In  Figure  6.7,  the  conversion  of  methanol  is  shown  as  a  function  of  methanol 
liquid  hourly  space  velocity  (LHSV-M)  for  conditions  with  and  without  an  acoustic  field. 
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Figure  6.6:  Difference  in  conversion  as  a  function  of  bed  length  for  4.8  and  9.5  g/min 

premix  flow  rates. 

In  this  figure  space  velocity  was  changed  by  adjusting  the  flow  rate  in  a  0.576  m  catalyst 
bed.     Note  that  there  is  a  significant  increase  in  the  conversion  for  those  cases  utilizing 
acoustic  waves.  This  figure  shows  that  one  can  increase  the  average  conversion  by  up  to 
4%  over  the  non-oscillating  cases  for  the  flow  rates  shown.  This  difference  in  conversion 
increases  for  higher  space  velocities  as  previously  shown  in  Table  6.1  and  Figure  6.6. 
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Figure  6.7:  Conversion  versus  space  velocity  with  and  without  a  165  dB  acoustic  field. 
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Notice  that  in  Figure  6.7  this  increase  in  conversion  by  using  an  acoustic  field  reduces  the 
wasted  methanol  by  approximately  50%  at  2.0  LHSV-M. 

At  a  given  conversion  the  space  velocity  can  be  increased  by  utilizing  the  acoustic 
field.  For  example,  in  Figure  6.7  at  a  96%  conversion  the  allowable  space  velocity 
increases  from  approximately  1.5  hr"1  to  2.0  hr"1  using  the  linear  regression.  This  is  an 
increase  of  33%  in  generating  capacity  or  reduction  in  reactor  volume  over  the  non- 
oscillating  case.  While  96%  fuel  conversion  still  may  be  an  intolerable  amount  of  fuel 
loss,  this  demonstrates  that  the  acoustic  field  can  increase  the  capacity  of  the  reformer. 
At  the  lower  end  of  space  velocities  tested  average  conversion  increases  of  0.6%  are 
actualized  by  employing  an  acoustic  wave.  While  this  does  not  seem  overly  significant  as 
a  conversion  increase,  the  amount  of  wasted  fuel  is  reduced  by  a  factor  of  2.5.  This 
increase  in  conversion  at  the  lower  end  of  space  velocities  becomes  more  important 
considering  that  certain  types  of  fuel  cells  can  be  poisoned  by  unconverted  methanol 
(Amphlett  et  al.,  1996).  While  the  conversion  increase  is  the  most  significant  output 
parameter  found  using  acoustics,  other  benefits  are  also  demonstrated. 

The  Effect  on  Controllability 

The  increased  heat  and  mass  transfer  due  to  the  acoustic  field  also  positively 
affects  the  controllability  of  temperature  within  the  reactor  for  all  locations.  Increased 
controllability  is  evidenced  by  both  a  reduction  in  the  variation  of  the  centerline 
temperature  and  an  decrease  in  the  average  difference  from  the  control  setpoint. 
Centerline  temperatures  were  used  as  these  locations  are  the  furthest  away  from  the 
feedback  thermocouple  positions  at  the  wall.  These  locations  showed  the  largest  average 
difference  from  the  setpoint  and  highest  variance  of  all  radial  locations  tested.  The  data 
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in  Table  6.2  show  average  and  standard  deviation  of  centerline  temperature  of  the  full 

length  catalyst  bed  at  two  premix  flow  rates  with  and  without  a  165dB(re20uPa)  acoustic 

field.  Further  data  demonstrating  enhanced  control  of  the  average  temperature  with 

Table  6.2:  Average  and  standard  deviation  of  the  centerline  temperatures  with  a  setpoint 
temperature  of  250°C  at  various  locations  both  with  and  without  a  165dB  acoustic  field. 


1.0LHSV-M 

2.0  LHSV-M 

With  Acoustics 

No  Acoustics 

With  Acoustics 

No  Acoustics 

Centerline 
Location 

Avg.  °C 

S.D. 

Avg.°C 

S.D. 

Avg.°C 

S.D. 

Avg.cC 

S.D. 

Zone  1 
Exit 

225.2 

2.12 

214.2 

1.87 

224.5 

1.41 

218.9 

1.67 

Zone  2 
Exit 

231.1 

2.70 

221.4 

3.12 

227.6 

1.38 

221.8 

2.14 

Zone  3 
Exit 

245.0 

3.10 

236.4 

4.13 

240.9 

1.99 

231.8 

2.42 

Zone  4 
Exit 

256.6 

4.87 

244.0 

2.96 

242.4 

1.89 

230.2 

1.48 

Zone  5 
Exit 

251.9 

1.95 

243.3 

1.04 

245.0 

2.12 

236.9 

1.52 

acoustics  was  presented  previously  in  the  temperature  profiles  as  Figure  6.1-6.4.  Clearly, 
an  acoustic  field  allows  for  better  control  of  temperature  throughout  the  reactor. 

The  Effect  on  Heat  Band  Power  Draw 
Related  to  the  change  in  controllability,  conversion  and  temperature  there  is  also  a 
change  in  the  power  demand  of  the  heat  bands.  Overall  average  power  draw  is  shown  in 
Table  6.3  and  the  use  of  an  acoustic  field  increased  the  overall  power  draw  slightly.  This 
accounts  for  an  overall  higher  temperature  and  increase  in  conversion  as  shown 
previously.  Although  the  overall  power  draw  was  similar  for  acoustic  and  non-oscillating 
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conditions  with  a  full  bed  length  (0.568  m)  at  two  space  velocities,  there  is  a  slight  shift 

in  where  the  power  is  drawn. 

Table  6.3:  Overall  average  catalyst  bed  power  draw  and  pooled  standard  deviation  for 
conditions  with  and  without  a  165  dB  acoustic  field. 


1.0LHSV-M 
No  Acoustics 

1.0LHSV-M 
With  Acoustics 

2.0  LHSV-M 
No  Acoustics 

2.0  LHSV-M 
With  Acoustics 

Pooled 
Standard 
Deviation 

158  Watts 

161  Watts 

212  Watts 

221  Watts 

Figure  6.8  shows  the  average  power  draw  in  each  zone  for  oscillating  (165  dB) 
and  non-oscillating  conditions  for  space  velocities  of  1.0  and  2.0  LHSV-M.  It  is  evident 
that  power  demand  at  all  zones  is  not  significantly  affected  by  the  acoustic  wave  but  the 
general  trend  is  that  the  power  draw  toward  the  entrance  of  the  oscillating  reactor  is 
higher  while  the  power  draw  closer  to  the  exit  is  higher  for  non-oscillating  conditions. 
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Figure  6.8:  Average  power  draw  for  zone  band  heaters  at  two  space  velocities  with  and 

without  a  165  dB  acoustic  field. 
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This  trend  indicates  that  the  acoustically  enhanced  reactor  is  undergoing  the  endothermic 
methanol  breakdown  further  upstream  than  in  the  non-oscillating  case.  Thus  for  the 
acoustically  enhanced  reformer,  there  is  more  space  and  time  for  the  fuel  to  complete  the 
conversion  to  final  products  rather  than  pass  through  the  catalyst  bed  unreacted. 

The  Effect  on  Power  Output 
As  shown  previously  the  input  heat  power  is  changed  slightly  by  employing  an 
acoustic  field.  Additionally,  the  acoustic  driver  also  requires  power.  Output  power  in 
terms  of  hydrogen  is  also  affected  by  the  acoustic  wave.  The  power  (VI)  required  by  the 
driver  for  various  output  sound  pressures  in  a  0.568  m  long  catalyst  bed  with  flowing 
reactants  is  shown  in  Figure  6.9.  For  a  pressure  amplitude  output  of  5300  Pa,  9.3  Watts 
(VI)  is  required.  The  power  factor  for  the  speaker  was  measured  to  be  0.92  at  650  Hz  so 
the  true  power  is  slightly  lower  than  the  9.3  Watts  measured.    The  output  power  in  Watts 
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igure  6.9:  Average  maximum  pressure  as  a  function  of  driver  power  at  650  Hz  with  a 
0.568  m  catalyst  bed  length  while  generating  hydrogen. 
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of  hydrogen  (LHV)  can  then  be  compared  to  this  input  power  to  determine  if  the 
enhancement  scheme  is  worthwhile.  As  premix  flow  rate  control  was  not  precise,  two 
representative  points  with  the  same  measured  premix  flow  rates  were  compared.  The 
total  average  dry  product  composition  was  assumed  to  be  representative  for  this 
comparison.  For  these  specific  experiment  runs  at  a  space  velocity  of  2.0  LHSV-M,  the 
dry  product  gas  flow  rate  was  found  to  increase  by  0.333  grams  per  minute  for  the  case 
with  acoustics.  With  the  average  hydrogen  mass  percentage  of  products  being  1 1 .81%, 
this  translates  into  an  increase  of  78.9  Watts  of  hydrogen  (LHV)  that  was  effectively 
enabled  by  9.3  Watts  (VI)  of  input  power  to  the  acoustic  driver  and  9  Watts  of  additional 
heat  demand.  This  yields  a  power  benefit-cost  ratio  of  over  4.  At  the  lower  flow  rates 
investigated  the  benefit-cost  ratio  decreased  to  about  1 .4.  Understanding  that  slight 
variation  in  measured  premix  flow  rates  did  exist,  Table  6.4  shows  the  averages  of 
additional  heat  demand,  the  input  power  of  the  driver,  the  dry  product  gas  flow  rate,  the 
hydrogen  (LHV)  output  power,  and  the  power  benefit-cost  ratio  for  different  space 
velocities  and  acoustic  conditions. 

It  is  important  to  note  that  the  power  required  was  the  power  (VI)  actually 
delivered  to  the  speaker  and  does  not  include  inefficiencies  due  to  the  generation  and 
amplification  of  the  reference  signal.  Also  note  that  the  reformer  has  not  been  optimized 
to  take  full  advantage  of  the  acoustic  wave  and  thus  greater  benefits  may  be  had  for 
acoustically  designed  reactors  or  monolith  catalysts  that  would  not  damp  out  the  acoustic 
wave.  Figure  6.10  shows  the  required  speaker  power  (VI)  for  attaining  sound  pressure 
amplitudes  in  the  0.095  m  length  catalyst  bed  or  one  sixth  of  the  available  catalyst  bed 
length.  As  shown,  it  was  found  that  the  driver  power  required  for  the  nominal  amplitude 
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Table  6.4:  A  summary  of  significant  performance  parameters  with  and  without  acoustic 

enhancement  for  two  space  velocities. 


Averages 

Additional 

Heat 
Demand 
(Watts) 

Acoustic 

Driver 

Power 

(VI) 

(Watts) 

Dry 
Product 

Flow 
(g/min) 

Hydrogen 
Power 
Output 
(LHV) 
(Watts) 

Benefit 
Cost 
Ratio 

2.0  LHSV-M  with  a 
165dB  Acoustic  Field 

9 

9.3 

7.58 

1792 

4.48 

2.0  LHSV-M  without 
an  Acoustic  Field 

N/A 

0 

7.23 

1710 

N/A 

1 .0  LHSV-M  with  a 
165dB  Acoustic  Field 

3 

9.3 

3.93 

929 

1.46 

1.0  LHSV-M  without 
an  Acoustic  Field 

N/A 

0 

3.85 

911 

N/A 

of  5300  Pa  could  be  obtained  with  much  less  power  when  the  catalyst  bed  length  was 
lowered.  This  is  due  to  decreased  acoustic  resistance  through  the  shortened  catalyst  bed 
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igure  6.10:  Average  maximum  pressure  as  a  function  of  driver  power  at  650  Hz  with  a 
0.095  m  catalyst  bed  length  while  generating  hydrogen. 
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length  and  is  evidence  that  acoustically  designed  catalyst  beds  may  have  much  higher 
benefit-cost  ratios  than  the  reactor  studied. 

The  Effect  on  Catalyst  Degradation 
Both  with  and  without  the  acoustic  field  the  catalyst  degrades  over  time.  This 
degradation  may  be  a  result  of  fouling  or  film  coating  of  the  catalyst  pellets  or  it  may  be 
due  to  oxidation  of  the  copper-zinc  catalyst  by  oxygen  present  in  the  reaction.  Another 
reason  for  degradation  may  be  that  of  sintering  of  the  catalyst  as  there  are  spatial  and 
temporal  gradients  through  the  reactor.  Acoustics  has  not  been  shown  to  significantly 
decrease  or  increase  the  degradation  rate  of  the  catalyst.  Figure  6.1 1  shows  the 
conversion  over  50  hrs  time  both  with  and  without  a  165  dB  acoustic  field  at  an  average 
premix  flow  rate  of  9.3  g/min  with  a  0.568  m  catalyst  bed  (-2.0  LHSV-M).  This  figure 
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Figure  6.1 1 :  Catalyst  degradation  with  and  with  out  acoustics  over  a  50  hour  time 
period  with  an  average  premix  flow  rate  of  9.3  g/min. 


99 

proves  that  utilizing  an  acoustic  field  in  the  reformation  process  has  little  or  no  effect  on 
catalyst  degradation  rate  for  the  conditions  studied.  However,  if  catalyst  life  can  be  seen 
as  time  or  the  amount  of  fuel  processed  above  a  set  limit  of  conversion,  acoustics  can 
help  the  process  because  the  conversion  starts  out  at  a  higher  level.  Utilizing  a  90% 
conversion  limit  for  the  operating  conditions  above.  Figure  6.1 1  shows  an  increase  of 
lifetime  or  the  amount  of  fuel  processed  by  2.5  times  for  cases  with  the  acoustic  field. 
Although  degradation  rate  was  not  observed  to  change  with  the  acoustic  field 
there  are  several  mechanisms  by  which  the  degradation  could  possibly  be  affected  by  the 
acoustic  waves.  It  has  been  shown  that  the  acoustic  field  decreases  the  spatial  and 
temporal  temperature  gradients  through  the  catalyst  bed.  This  should  allow  for  decreased 
sintering  of  the  catalyst.  Furthermore  the  decreased  temperature  gradient  associated  with 
acoustics  would  not  allow  for  temperatures  approaching  the  dewpoint  of  the  mixture  to  be 
reached.  This  may  prevent  fouling  due  to  condensed  species  which  fouling  has  been 
shown  to  have  a  permanent  detrimental  effect  on  catalyst  performance  (Amphlett  et  al., 
2001 )  As  previously  discussed  in  the  theoretical  approach  (Chapter  2).  the  recirculation 
of  hydrogen  due  to  acoustic  mixing  may  also  alleviate  poisoning  through  oxidation  of 
reduced  catalyst.  In  order  to  determine  the  effect  of  poisoning  due  to  possible  oxidation 
of  the  catalyst,  data  showing  recovery  of  the  catalyst  was  gathered.  After  the  periods  of 
operation  shown  in  Figure  6.1 1,  the  catalyst  bed  was  subjected  to  a  volume  fraction  of 
2.3%  hydrogen,  97.7%  nitrogen  mixture  for  a  period  of  1 0  hours  at  1 80 °C.  This  was 
followed  by  an  hour  at  220"C  increasing  the  volume  fraction  of  hydrogen  to  6.5  %. 
Another  hour  at  220  °C  temperature  with  a  12%  volume  of  hydrogen  followed.  This 
effectively  re-reduced  the  catalyst  and  would  reverse  any  poisoning  due  to  oxidation. 
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However,  recovery  of  the  methanol  conversion  was  not  observed  for  either  case  with  or 
without  the  acoustic  field.  This  indicates  that  the  degradation  of  the  catalyst  bed  is  due  to 
fouling,  sintering  or  poisoning  other  than  simple  oxidation. 

The  degradation  demonstrated  is  not  necessarily  typical  of  methanol  steam 
reforming  in  general  but  is  specific  to  the  conditions  tested.  The  degradation  studies  took 
place  at  double  the  original  design  flow  rate  of  fuel  to  show  aging  of  the  catalyst  within  a 
reasonable  period  of  time.  Additionally  impurities  may  have  come  into  the  catalyst  bed 
with  the  industrial  grade  methanol  fuel.  Whatever  the  specific  degradation  mechanisms 
in  the  studied  reactor,  employing  an  acoustic  field  has  been  shown  to  have  little  or  no 
effect  on  catalyst  degradation  rate. 

Potential  Benefits  for  a  Fuel  Cell  Vehicle 

The  increase  in  conversion  for  a  given  space  velocity,  and  correlated  increased 
power  output  translate  into  a  direct  benefit  for  fuel  cell  vehicle  applications.  Utilizing 
acoustic  fields  may  allow  a  feasible  method  of  decreasing  the  volume  and  weight  of  the 
reformer  through  the  increase  in  reformer  capacity  and  shown  in  Figure  6.7  and  Table  6.1. 
Transients  at  start-up  would  also  be  influenced  by  the  ability  to  have  a  smaller  reactor 
even  without  including  the  influence  of  higher  heat  transfer  within  the  reactor  due  to 
acoustics.  This  is  shown  by  considering  a  simplified  energy  balance  around  the  reactor 
for  a  start-up  scenario.  Simplifying  the  process  by  assuming  that  heat  will  be  transferred 
into  the  catalyst  bed  by  convection  only  and  that  the  reactor  will  have  minimal 
temperature  gradients  throughout  (lumped  capacitance  model),  an  energy  balance  yields 
the  differential  equation  in  Equation  6. 1 .  This  is  a  model  for  the  temperature  response  of 
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hAT+T=T"f  (6'1} 

the  reactor  subjected  to  a  step  function  in  temperature  (from  T0  to  T,nt).  as  would  be  a 

simplified  case  with  start  up.  Equation  6.2  shows  the  solution  of  this  ordinary  differential 

equation  and  yields  the  time  response.  It  is  clear  that  as  one  decreases  the  mass  of  the 

reactor  the  response  time  is  decreased  proportionally.  The  effect  could  be  further 

magnified  by  considering  the  increased  convective  heat  transfer  rate  within  the  reformer. 

T(t)  =  (T0-Tinf)exp  -—   +Tinf  (6.2) 

Considering  the  data  presented,  it  may  seem  that  acoustics  is  a  clear  choice  for 
improving  the  performance  of  a  steam-reforming  reactor  in  a  fuel  cell  vehicle.  However, 
there  are  drawbacks  associated  with  practically  superimposing  the  sound  waves  in  the 
steam-reforming  processor.  With  the  current  system  the  acoustic  wave  must  be  generated 
by  an  acoustic  driver.  If  implemented  in  a  fuel  cell  vehicle,  this  driver,  along  with  an 
amplification  system  and  control  would  add  to  the  volume,  weight  and  complexity  of  the 
vehicle.  For  the  system  studied,  the  relationship  between  auxiliary  power  draw  (Paux)  and 
electrical  power  input  to  the  driver  (Pdriver)  at  650  Hz  was  found  as  Equation  6.3.  For  the 
current  system  the  auxiliary  power  draw  to  generate  the  acoustic  waves,  along  with  the 

PA,„  =  2.53(PJmJ+ 117.4     (Watts)  (6.3) 

weight  and  volume  of  the  driver  is  shown  in  Table  6.5.  The  potential  benefits  of  the 
acoustical  enhancement  must  be  balanced  against  these  practical  drawbacks  before 
implementation  in  a  fuel  cell  vehicle  is  justified. 
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Table  6.5:  Auxiliary  requirements  for  generating  the  acoustic  fields  in  the  current  facility. 


Auxiliary  Power  Requirement 

(Amplification  and  Signal 

Generation) 


Driver  and 
Housing  Weight 


Driver  and  Housing 
Volume 


140  W  (9  W  to  driver) 


3.9  kg 


6630  cm3 


Potential  decreases  in  time  response  at  start-up  due  to  the  acoustic  fields  also 
should  be  viewed  in  light  of  the  rate  limiting  step  in  the  transient  response.  During 
start-up  of  an  fuel  cell  bus  as  described  by  Erickson  et  al.  (2001),  the  initial  heating  of  the 
reformer  is  limited  by  catalyst  overheating  at  the  walls  due  to  the  temperature  gradient 
inside  of  the  reactor.  This  temperature  gradient  may  be  changed  significantly  by  the 
acoustic  wave  as  described  previously,  but  most  catalysts  typically  have  some  type  of 
limit  in  heating  rate  to  avoid  sintering.  If  this  catalyst  heating  rate  is  less  than  the 
possible  heating  rate  of  the  overall  reformer  when  considering  start-up,  a  smaller  and 
faster  heating  reformer  due  to  acoustic  waves  will  not  assist  with  this  type  of  transient 
operation.  However,  transient  operation  that  does  not  require  a  change  in  overall 
temperature  such  as  with  load  following  may  be  enhanced  by  the  acoustic  wave  due  to  the 
increase  in  the  heat  transfer  rate  within  the  reactor  as  presented  previously  in  Chapter  2. 
Further  data  is  needed  to  quantify  the  enhancement  and  limitations  of  using  acoustic 
waves  for  transient  operation. 

Although  the  conversion  is  increased  by  employing  acoustic  fields,  and  the  power 
benefit-cost  ratio  is  favorable  when  considering  power  into  the  reformer,  there  are 
systems  that  need  to  be  optimized  and  integrated  before  practical  implementation  is 
feasible.  The  primary  concerns  are  the  size  and  weight  of  the  acoustic  driver  system,  and 
power  demand  of  the  reference  signal  generation  and  amplification. 


CHAPTER  7 
CONCLUSIONS  AND  RECOMMENDATIONS 

From  the  theory  and  data  presented  previously,  conclusions  can  be  drawn  and 
recommendations  for  further  study  and  development  can  be  made 

Conclusions 

A  theory  was  developed  that  points  toward  enhancement  of  the  steam-reforming 
process  by  utilizing  an  acoustic  field  in  critical  fluid  flow  paths.  Specific  mechanisms 
and  models  of  this  enhancement  were  reviewed.  An  acoustic  field  is  thought  to  enhance 
the  reformation  process  primarily  through  increasing  the  heat  transfer  and  mass  transfer 
properties  within  the  catalyst  bed.  An  acoustic  field  generator  was  attached  to  a  methanol 
steam  processor  and  utilized  to  control  the  acoustic  field  within  the  catalyst  bed.  The 
sound  pressure  level  made  possible  by  this  generator  is  165  dB  (reu20Pa)  at  650  Hz  while 
producing  hydrogen.  Empirical  transfer  functions  indicate  other  available  resonant 
frequencies  but  at  a  lower  sound  pressure  levels.  The  data  collected  on  conversion  of 
methanol  has  been  worked  into  an  empirical  model  that  accurately  predicts  the 
conversion  for  a  set  of  acoustic  pressure  amplitude,  and  space  velocity  parameters.  The 
acoustic  field  affects  several  measurable  parameters  within  the  catalyst  bed.  As  a  result  of 
the  increased  heat  and  mass  transfer  found  with  acoustics,  the  radial  temperature  gradient 
within  the  catalyst  bed  becomes  more  even  but  is  not  eliminated.  Similarly,  the  reactor 
was  shown  to  be  better  controlled  with  an  acoustic  field  in  average  centcrline  temperature 
and  in  deviation  from  the  average  for  most  axial  positions.  The  conversion  of  methanol 
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fuel  was  increased  significantly  by  employing  an  acoustic  field.     This  increase  in 
conversion  allows  a  reduction  of  required  catalyst  bed  length  or  an  increase  in  capacity 
for  a  given  reformer  volume  with  an  acoustic  field.  Empirical  data  show  that  a  25% 
decrease  in  reformer  volume  or  a  similar  increase  in  flow  capacity  can  take  place  by 
employing  a  resonant  acoustic  field  within  the  methanol-steam  reformer.  The  power 
output  in  hydrogen  was  likewise  increased  by  utilizing  an  acoustic  field.  The  benefit-cost 
ratio  of  using  acoustics  in  terms  of  power  ranges  from  1.7  to  4  for  the  flow  rates  and  bed 
lengths  studied.  This  favorable  benefit-cost  ratio  points  toward  practical  implementation 
of  acoustically  enhanced  reformers  in  applications  such  as  fuel  cell  vehicles.  Catalyst 
degradation  was  not  notably  affected  by  use  of  an  acoustic  field.  All  of  these 
characteristics  showing  enhancement  of  the  process  were  found  without  optimizing  the 
catalyst  bed  or  acoustic  generator  for  use  in  this  application.  These  results  prove  that 
acoustic  enhancement  of  reformation  process  can  occur  under  the  proper  conditions. 

Recommendations 
Although  this  study  has  demonstrated  that  use  of  acoustics  is  a  viable 
enhancement  method,  there  are  further  recommended  studies  and  applications.  The 
acoustic  driver  system  was  not  optimized  for  this  application  but  was  constructed  from 
"off  the  shelf  available  parts.  Likewise,  the  steam  reformer  was  not  designed  to 
maximize  the  propagation  of  the  acoustic  waves  and  the  pelletized  catalyst  was  shown  to 
significantly  limit  the  acoustical  resonance  within  the  reactor.  Considering  these 
limitations,  this  study  proves  that  reforming  in  conjunction  with  acoustics  is  feasible  and 
significantly  influences  the  reformation  process  even  under  non-ideal  circumstances. 
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Understanding  that  such  feasibility  exists,  reformers  and  acoustic  drivers  should  be 
constructed  to  optimize  the  proven  enhancement  by  using  acoustic  fields. 

The  capability  of  higher  sound  pressure  levels  within  the  reactor  has  been  proven 
and  an  increase  in  sound  pressure  level  should  increase  the  observed  benefits,  further  than 
shown  in  this  study.  Higher  sound  pressure  levels  may  be  possible  by  utilizing 
monolithic  catalysts,  reducing  the  acoustic  path  in  the  reactor,  or  utilizing  a  different 
actuating  system. 

As  presented  in  Chapter  2.  monolithic  catalysts  should  be  affected  by  the  acoustic 
wave  substantially  more  than  the  pelletized  catalyst  used,  even  without  any  sound 
pressure  level  increase.  This  is  due  to  the  dominance  of  internal  area  with  pelletized 
catalysts  which  is  not  considered  to  be  affected  by  the  periodic  changes  in  gas 
displacement  at  the  surface  associated  with  the  acoustic  waves.  Data  utilizing  monolithic 
catalysts  should  be  gathered  and  compared  to  the  enhancement  found  in  this  study. 

Reformers  and  acoustic  drivers  should  also  be  optimized  for  utilizing  the  velocity 
anti-nodes  of  the  acoustic  modes.  This  could  include  the  use  of  resonant  chambers, 
flapper  valves,  acoustic  filters,  or  removal  of  the  catalyst  not  affected  by  the  velocity  anti- 
nodes.  Additionally  a  movable  catalyst  bed  within  an  acoustically  enhanced  reformer 
may  be  able  to  yield  the  optimal  sections  for  enhancement  that  were  alluded  to  in  this 
study. 

Ideal  reactor  geometry  should  also  be  investigated  as  to  temperature.  This  stems 
from  the  fact  that  temperature  gradients  have  been  shown  to  be  reduced  by  the  acoustic 
field  but  not  eliminated.    Combinations  of  acoustics  and  more  conventional  mixing 
enhancement  such  as  induced  swirl  should  be  investigated.  Furthermore,  the  reforming 
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enhancement  from  acoustic  waves  should  be  compared  and  contrasted  to  conventional 
methods. 

As  acoustic  driver  mass  and  volume  and  overall  power  draw  need  to  be 
minimized,  system  optimization  should  occur.  An  overall  acoustic  system  needs  to  be 
designed  for  this  specific  application  because  the  power  draw  of  the  amplifier  and 
generator  system  is  currently  too  high  to  justify  implementation.  As  endothermic  and 
exothermic  regions  exist  within  the  catalyst  bed,  natural  Rijke-type  oscillations  are 
theoretically  possible  and  thus  may  eliminate  the  need  for  auxiliary  components  and 
drivers.  This  possibility  should  also  be  investigated. 

Transient  operation  of  the  reformer  with  acoustics  should  be  studied.    Further 
data  is  needed  to  quantify  the  potential  benefits  of  decreased  time  response  during 
transient  operation.  This  includes  transients  at  start-up  and  load  following. 

Autothermal  and  partial  oxidation  reforming  are  other  hydrogen  production 
methods  that  may  be  influenced  by  acoustic  waves.  This  would  be  particularly  interesting 
in  autothermal  methods  because  the  exothermic  portions  of  the  catalyst  bed  could 
combine  with  the  partial  recirculation  effects  of  resonant  acoustic  waves. 

In  summary,  this  study  has  proven  that  the  acoustic  waves  are  a  feasible  method  of 
enhancing  the  reformation  of  hydrocarbon  fuels  into  hydrogen.  Further  studies,  including 
those  recommended  above,  may  enable  practical  implementation  of  this  technology  for 
fuel  cell  vehicles. 


APPENDIX  A 
COMPUTER  PROGRAMS 

Several  computer  programs  were  used  to  collect  data  from  the  instruments. 

control  the  temperature,  monitor  the  acoustic  field,  and  analyze  the  data.  These  programs 

were  built  using  the  Graphical  User  Interface  (GUI)  of  National  Instruments'  LabVIEW. 

The  LabVIEW  Programs  used  in  this  study  are  displayed  in  the  following  pages.  These 

include  the  front  panel  (user  interface)  and  block  diagram  (logic)  for  each  of  the  programs 

used.  The  first  program  presented  is  the  temperature  control  and  data  acquisition 

program  utilizing  IOTech  DBK-52  and  DBK-24  hardware  named  Controller.  Following 

this  is  the  program  called  Spectrum. vi  which  utilizes  a  National  Instruments  DAQCard- 

1 200  to  collect  the  data  and  produced  the  transfer  functions  presented.  Next  the  data 

acquistion  program  named  PaulsScope  is  shown.  This  program,  using  the  DAQ-Card 

1 200  analyzed  the  raw  voltage  traces  from  the  pressure  transducer,  translated  it  into 

maximum  sound  pressure  and  sound  pressure  level,  and  saved  these  data  to  disk. 

Analysis  programs  follow  these  data  aquisition  programs.  The  first  such  program  is 

called  Conversion.  This  program  utilized  the  amount  of  premix  used,  amount  of 

condensate  collected  and  condensate  density  to  calculate  the  methanol  conversion.  The 

next  program  called  Space  Velocity  calculated  both  the  LHSV-M  and  the  Gaseous  Hourly 

Space  Velocity  (GHSV)  for  a  set  of  temperature  flow  rate  and  catalyst  bed  parameters. 
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Controller  Front  Panel 


Controller  Block  Diagram 


oonapDappoooooaaoooopn  out  n  n ,  b]k|Dppppnoppppppppppppppop 


500  00 


B7E 


|Tempeialuiei|   jjg), 

(Time    inr  r  l.i  I  [■■  --ij 


a 


Irrril 


shuttle  Inpuls  lo  conecl  output  form 
land  check  surfaces 


E> 


E> 


sm     ■* 


E> 


T'ar 


'"'V     -if--.!- 


..£a> G 


jQdqti 


□  y 


HBHoannapgn  nTnfrrp  mm  otto  flBBBBBBgngflBaBaBBBBannnggaoo-o' 


Controller  Block  Diagram  Continued 


109 


QOonnnnnoongnaoDnnDaD fl_q<|  i  [u  t:i^n DnnaDonoagaoannnannptn 

(BiddContiolondTA*! 


@HPF~  ~mji 


E> 


tfr 


E> 


E> 


IthUfe  Inputs  to  coned  output  lorm| 

Idtsptay  lempetatmei  anddT/c>  contipj   GfTS] 


T^P       f 

Vapors 

HE  t 

-• 

■ 

MB  4 

ne  5 

HS  ?. 

i;...-:l^.-|.- 

i> 


nF3 


3    Kataialt  itopelj 


CD 


t> 


add  n  n 


1!R>^ 


q-n-nonnnnnnnnnnBnnnntiBnBBnnHnanngBmt 


anpoDOODPODogpoDPoan OJum  2|0  6iHDflpaoHogpQqponaoooo.gotn 


""Hl^E" 


aiti 


BpJ 


j  fm  On]  onbol 


MH>- 


CD 


Tkiw  ( in  i  I'JnJ) 


tunnm: 


n  n  b  n  a  a  mnra-nrmnrt 


1*5' 


1QE) 


EQ 


Controller  Block  Diagram  Continued 
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Spectrum  Block  diagram 
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Spectrum  Block  diagram  continued 
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Spectrum  Block  diagram  continued 
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Spectrum  Block  diagram  continued 
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PaulsScope  Front  Panel 
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PaulsScope  Block  Diagram  Continued 
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Space  velocity  Front  Panel 
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APPENDIX  B 
CHECKLISTS  AND  PROCEDURES 

The  following  pages  are  the  checklists  used  for  start-up.  run  (or  data  collection). 

and  shutdown.  These  checklists  were  generated  using  previous  studies  performed  on  the 

reactor  and  experience  gained  through  preliminary  investigation. 
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24 


Start-up  Date Time 

□  Prepare  Water-Methanol  Premix  for  1 .5: 1  adjust  to  .9075-.9085 

(A)  Empty (B)MeOH+Empty (B-A)= *. 84375= 

+B= 

□  Setup  Tubing  in  Pump.  Pump  Setpoint ml/min 

□  Switch  vaporizers  to  on  at  power  box 
Check  valves  E,&  D  open,  CO,  open,  cooling  water  on,  set  N2  flowrate,  R,&  C  closed. 


T 

dT/dt 

Max  time 

Vaporizer  Surface 

420 

1.5 

90 

Vaporizer  Exit 

250 

99 

Trap  Surface 

480 

1.5 

90 

Trap  Exit 

250 

99 

Superheater  Surface 

275 

1.5 

80 

Superheater  Exit 

250 

.5 

Zone  1  -5  Exit 

250 

0.3 

40 

Zone  6  Exit 

170 

0.3 

40 

□  Switch  Relay  Timer  to  1 00  (Always  On),  RESET 

□  Turn  on  Vaporizers  and  after  Vaporizer  surfaces  reach  1 50 °C  turn  on  Pump  at  Control 
Panel 

□  Start  GC, 

□  Carrier  On 

□  Power  On  and  Check  Carrier  pressure 

□  TCD  High  for  He,  Low  for  N2 

□  H2  on  and  Check 

□  Air  on 

□  Ignite  FID  and  Check 

□  Run  once  with  Air 
When  trap  exit  reaches  250 °C 


T 

dT/dt 

Max  time 

Vaporizer  Surface 

420 

1.5 

Off 

Vaporizer  Exit 

280 

99 

Trap  Surface 

480 

1.5 

Trap  Exit 

250 

99 

Superheater  Surface 

305 

1.5 

Superheater  Exit 

250 

.5 

Zone  1-5  Exit 

250 

0.0 

Zone  6  Exit 

250 

0.0 

□  Open  D,  close  CO,,  open  R,  close  E 

□  Set  Acoustic  Pressure  and  Frequency  and  adjust  for  resonance 
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RUN  Date Time 

Catalyst  Weight 

□  Premix  Density g/ml  Water  Flow ml/min  N2  Flow 

□  Prepare  ice  bath  hook  up  condenser  and  turn  on 

□  (at  50  Min  or  stable  temperatures)  check  temp  (condenser  Should  be  0  deg) 

□  Write  Acoustic  Pressure (A volts/2)  /  1 05.6e-3= Pa  (Amplitude) Hz 


Acoustic 


filename 


□  Change  write  Interval  to sec,  write  file  on 

□  RESET  Open  C,  Close  D,  read  (1)  initial  premix  mass 

□  Start  Run  time  Filename 


□  Collect  GC  samples  at , ,  and , . 

□  H2,CO  and  C02  levels  H2 CO CO, H2 CO CO,_ 

□  At  end  of  run  off  write  file  Open  D,  Close  C,  read  (F)  final  premix  mass 

□  F-I= 

□  End  Run  time 

□  Collect  Condensate,  Weigh  (Z)condensate  +  empty g 

□  Weigh  empty  condenser(X) g    Z-X= g 

□  Power  measurements  Driver V (if  available)  Amps 


ANALYSIS 

□  Condensate  Density  at  24C g/ml 

□  Premix  Flowrate g/min ml/min 

D  Space  velocity _MeOH  GHSV 

□  Conversion 
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SHUTDOWN 

□  Ramp  Set  point  to  3  ml/min  (neglect  if  lower) 

□  Turn  off  condenser  and  drain 

□  Turn  off  Write  to  File 

at  10  min  after  ramp  down 

D  Open  E  Close  R  open  CO,  Open  C  Close  D 

□  Switch  off  vaporizer  and  trap  at  power  box 

□  Turn  off  Vaporizer  Trap  and  Pump  at  control  panel 

□  Switch  off  Pump  unlatch  tubing 

□  Relay  timer  to  1-10 

□  Switch  off  Breaker 

□  Switch  on  Breaker  (if  continuing  to  hold  at  temperature)  and  check  relay  timer 


T 

dT/dt 

Max  time 

Vaporizer  Surface 

NA 

NA 

off 

Vaporizer  Exit 

NA 

NA 

Trap  Surface 

NA 

NA 

Trap  Exit 

NA 

NA 

Superheater  Surface 

250 

10 

Superheater  Exit 

250 

10 

Zone  1-5  Exit 

250 

10 

Zone  6  Exit 

170 

10 

□  set  N2  flow  rate  to  30 

□  when  temperatures  reach  255  °C  turn  off  C02 

□  Shutdown  GC 

□  exit  program 

□  TCD  current  off 

□  Air  off 

□  H2  off 

□  Power  off 

□  Carrier  off 


APPENDIX  C 
STATISTICAL  DESIGN  OF  EXPERIMENTS 

This  appendix  is  designed  to  review  the  use  of  statistical  design  of  experiments  as 
utilized  in  Chapter  5  to  develop  a  statistically  valid  model.  This  is  the  authors 
simplification  and  summary  of  work  presented  by  Lawson  and  Erjavec  (2000). 

Factorial  design  of  experiments  is  a  method  of  building  an  empirical  model 
covering  a  design  space,  or  variable  space  using  the  least  amount  of  data  points.  Each 
input,  or  independent  variable  is  considered  at  endpoints  of  low  value  and  a  high  value 
corresponding  to  -1  and  +1,  respectively.  The  output  or  dependent  variable  is  measured 
at  the  corresponding  endpoints  of  the  design  space.  The  use  of  the  endpoints  allows  one 
to  determine  the  required  number  of  experimental  trials.  One  calculates  the  minimum 
required  number  of  trials  (Trialsmin)  by  raising  the  levels  (number  of  endpoints)  (V)  to  the 
power  of  the  number  of  independent  variables  (k)  as  shown  in  Equation  C.l .  For 

Trialsmin  =  Vk  (C.l) 

example  if  2  independent  variables  were  studied  and  each  variable  had  2  levels  (2 
endpoints)  the  minimum  required  number  of  trials  would  be  4.  This  is  called  a  22 
factorial  experiment.  The  resulting  factorial  matrix  would  appear  as  Tabic  C.l. 
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Table  C.l :  Factorial  matrix  for  2  input  variables  tested  at  2  levels  (+1  and  -1 ). 

Run 

Input  Variable  1 

Input  Variable  2 

Output  Variable 

1 

+  1 

+  1 

Y.  , 

2 

+  1 

-1 

Y+. 

3 

-1 

+  1 

Y  , 

4 

-1 

-1 

Y  . 

One  can  note  that  the  design  space  for  the  levels  of  the  input  variables  is  covered 
by  the  output  variables  collected,  or  in  other  words  the  output  is  known  for  all  individual 
variables  or  combinations  of  the  two  input  variables.  The  change  that  a  variable  makes 
on  the  output  is  known  as  the  effect  of  that  variable.  The  change  that  a  combination  of 
two  or  more  variables  makes  on  the  output  is  known  as  the  interaction  of  those  variables. 
Assuming  a  linear  relationship  across  the  design  space  one  can  predict  the  output  within 
the  design  space.  However,  the  linear  relationship  should  be  verified  by  collecting  data 
near  the  midpoints  for  the  input  variables. 

This  prediction  of  how  an  output  will  respond  due  to  the  variables  investigated 
requires  the  construction  of  a  numerical  model  from  the  data  collected.  The  effects  are 
calculated  by  taking  the  average  output  when  a  variable  is  at  its  high  level  minus  the 
average  output  when  the  variable  is  at  its  low  level.  This  is  the  effect  of  that  variable. 
Interactions  are  found  in  a  similar  manner  by  taking  the  average  output  when  the 
combination  of  two  or  more  variables  is  high  minus  the  average  output  when  the  variable 
combination  is  low  (see  Equation  5.2).  These  effects  and  interactions  can  be  used 
together  with  the  input  variable  levels  and  average  output  to  predict  the  output  for  any 
combination  of  variables  within  the  design  space.  An  example  of  this  model  is  shown  in 
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Equation  C.2  for  a  22  factorial  experiment.  The  difference  between  the  predicted  value 
and  the  mean  values  for  any  set  of  variables  is  known  as  the  residuals  and  these  should 
not  vary  by  more  than  the  standard  error.  Individual  data  points  should  not  vary  from 
their  mean  value  by  more  than  twice  the  pooled  standard  deviation. 


^fx2++%x1x2 


-x7.il 

2 


E, 
2 


(C.2) 


The  run  order  shown  for  the  example  in  Table  C.l  is  called  the  standard  order  or 
Yates  order.  Randomization  of  the  run  order  is  recommended  to  account  for  process 
drift.  The  reason  for  this  is  that  process  drift  can  misrepresent  the  effects  that  an 
individual  variable  or  variables  have  on  an  output.  For  example,  if  one  considers  a  2? 
factorial  matrix  (which  requires  a  minimum  of  8  trials)  run  in  the  Yates  order  for  the  first 
input  variable  (X,)  with  process  drift  shown  by  Figure  C.l,  the  process  drift  misrepresents 


Q. 


Significant 
difference 


Y  + 


Time 


Figure  C.l :  Improperly  designed  experiment  without  a  randomized  run  order, 
the  change  in  output  due  to  X,  significantly  as  there  is  a  large  difference  between  average 

outputs  for  positive  values  of  X,  (Y+)  and  negative  values  of  X,  (Y-).  In  Figure  C.2, 
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where  the  run  order  has  been  randomized,  no  misrepresentation  occurs  because  there  is 
little  difference  between  the  average  output  when  the  variable  is  positive  (Y+)  and  the 
average  of  the  output  when  the  variable  is  negative  (Y-)  or  at  its  low  level. 


3 
Q. 

3 

o 


Time 


Difference  not 
significant 


Y  + 


Y- 


Figure  C.2:  A  properly  designed  experiment  with  a  randomized  run  order. 


Statistical  analysis  can  take  place  if  replications  are  made  of  the  study.  Each 
replication  increases  the  required  number  of  trials.  Thus  for  a  23  factorial  experiment 
with  2  replications,  24  individual  data  points  would  be  required.  These  replications  allow 
statistical  analysis  or  uncertainty  to  be  calculated  on  each  variable  setting  and  a  combined 
or  pooled  uncertainty  can  be  determined  for  the  data  set  and  for  the  resulting  model.  The 
uncertainty  for  each  setting  of  variables  can  be  calculated  from  the  resulting  standard 
deviation  and  a  t-table.  Overall  uncertainty  is  estimated  by  doubling  the  average  standard 
deviation  and  dividing  by  the  square  of  the  number  of  data  points  gathered  as  shown  in 
Equation  5.1.  This  uncertainty  estimate  is  known  as  the  standard  error  or  the  expected 
error  of  the  mean  value. 
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Because  experimental  errors  and  process  drift  are  always  present,  the  effects  and 
interactions  found  from  construction  of  the  model  must  be  tested  to  find  their  level  of 
significance.  If  a  high  level  of  noise  is  present  one  may  make  an  unwarranted  conclusion 
that  a  variable  has  a  large  effect  on  the  output  variable  when  the  variation  in  the  output  is 
only  due  to  the  noise  of  the  system.  The  calculation  of  the  signal-to-noise  t-ratio  allows 
comparison  to  standard  noise  levels  at  various  confidence  intervals.  The  signal-to-noise 
t-ratio  is  calculated  by  dividing  the  effect  or  interaction  by  the  standard  error  (see 
Equation  5.4).  Comparison  of  the  absolute  value  of  the  signal-to-noise  t-ratio  to  the 
critical  t-value  at  the  total  degrees  of  freedom  and  desired  confidence  level  allows  one  to 
determine  the  statistical  significance  of  the  effect  or  interaction.  Those  absolute  values  of 
signal-to-noise  t-ratios  greater  than  the  t-value  are  considered  significant  at  the 
confidence  level  investigated.  Those  which  are  smaller  are  not  considered  significant  at 
that  confidence  level. 

The  following  pages  are  an  example  of  a  23  factorial  experiment  with  one 
replication  to  investigate  the  accuracy  of  a  throwing  device.  This  is  adapted  from  a 
similar  example  given  by  Lawson  and  Erjavec,  (2000).  This  example  includes 
construction  of  the  model  and  a  statistical  test  of  the  signal  to  noise  t-ratios. 
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Table  C.2:  Hypothetical  Problem 

Factors  and  Levels 

Factor 

Levels 

+ 

- 

Xj=  Ambient  temperature 

32°C 

22°C 

X2=  Throwing  Machine  warmup  time 

5  min 

30  sec 

X3=  Time  that  power  is  connected 

5  min 

30  sec 

Y=  Measured  Distance  (m) 

Tab 

le  C.3:  Data  from 

Hypothetical  Problem 

Design  and  Data 

Run 

Coded  Values 

Actual  Values 

Run 
Order 

Y(m) 

XI 

X2 

X3 

x, 

X2 

x3 

1 

- 

- 

- 

22°C 

30  sec 

30  sec 

5,4 

69.2 

69.3 

2 

+ 

- 

- 

32°C 

30  sec 

30  sec 

14,7 

63.6 

63.5 

3 

- 

+ 

- 

22°C 

5  min 

30  sec 

15,2 

69.1 

69.4 

4 

+ 

+ 

- 

32°C 

5  min 

30  sec 

1,3 

63.1 

63.3 

5 

- 

- 

+ 

22°C 

30  sec 

5  min 

8,11 

69.1 

63.5 

6 

+ 

- 

+ 

32°C 

30  sec 

5  min 

12,16 

68.0 

67.9 

7 

- 

+ 

+ 

22°C 

5  min 

5  min 

10,6 

69.4 

69.3 

8 

+ 

+ 

+ 

32°C 

5  min 

5  min 

9,13 

68.4 

63.6 
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Table  C.4:  Calculation  c 

f  Averagt 

3S,  Variance  and  Degrees  of  Freedom 

Run 

x, 

x2 

x3 

Y 

Y 

avg 

Variance 

Degrees 

of 
Freedom 

1 

- 

- 

- 

69.2 

69.3 

69.25 

0.005 

2 

+ 

- 

- 

63.6 

63.5 

63.55 

0.005 

3 

- 

+ 

- 

69.1 

69.4 

69.25 

0.045 

4 

+ 

+ 

- 

63.1 

63.3 

63.20 

0.02 

5 

- 

- 

+ 

69.1 

63.5 

66.30 

15.68 

6 

+ 

- 

+ 

68.0 

67.9 

67.95 

0.005 

7 

- 

+ 

+ 

69.4 

69.3 

69.35 

0.005 

8 

+ 

+ 

+ 

68.4 

63.6 

66.00 

1 1 .52 

Table  C.5:  Pooled  Standard  Deviation  Total  Degrees  of  Freedom  and  Total  Average 


Pooled  Standard 
Deviation 

Total  Degrees  of 
Freedom 

Average  of  all  Ya 

1.847 

8 

66.86 

Table  C.6:  Worksheet  for  calculation  of  Effects  and  Interactions 


Run 

x, 

x2 

x3 

x,x2 

x,x3 

X2X3 

x,x2x3 

Yav8 

1 

- 

- 

- 

+ 

+ 

+ 

- 

69.25 

2 

+ 

- 

- 

- 

- 

+ 

+ 

63.55 

3 

- 

+ 

- 

- 

+ 

- 

+ 

69.25 

4 

+ 

+ 

- 

+ 

- 

- 

- 

63.20 

5 

- 

- 

+ 

+ 

- 

- 

+ 

66.30 

6 

+ 

- 

+ 

- 

+ 

- 

- 

67.95 

7 

- 

+ 

+ 

- 

- 

+ 

- 

69.35 

8 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

66.00 

From  pooled  standard  deviation  ( 1 .847)  and  total  number  of  data  points  (16)  Standard 
Error  is  calculated  to  be  0.9239. 
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Tab 

e  C.7:  Calculation  of  effects  and  signal-to-noise  t-ratios. 

x, 

X2 

x3 

X,X, 

XjX3 

X2X3 

X,X2X3 

E(+) 

260.7 

267.8 

269.6 

264.75 

272.45 

268.15 

265.1 

I(-) 

274.15 

267.05 

265.25 

270.1 

262.4 

266.7 

269.75 

£(+)-£(-) 

-13.45 

0.75 

4.35 

-5.35 

10.05 

1.45 

-4.65 

(Z(+)-I(-))/4 

or  Effects 

-3.36 

0.19 

1.09 

-1.34 

2.51 

0.36 

-1.16 

Signal-to-noise 
t-ratio 

-3.64 

0.20 

1.18 

-1.45 

2.72 

0.39 

-1.26 

Critical  t-ratio  at  95%  Confidence  and  8  degrees  of  Freedom  =  2.306 


Table  C.8:  Summary  of  statistical  significance  for  the  hypothetical  problem 

E, 

E2 

E3 

1,2 

1,3 

23 

1,23 

Statistically  Significant  ? 

Yes 

No 

No 

No 

Yes 

No 

No 

The  statistically  significant  model  of  the  design  space  is  predicted  by  Equation 
C.3  and  with  numerical  values  in  Equation  C.4 


Y=Y+|xi+^X1X3 

-3.36         251 
Y  =66.86  +  ^r— X1+^—X1X3 


(C.3) 
(C.3) 


APPENDIX  D 
RUN  ORDERS 

The  following  pages  are  the  set-ups  and  run  orders  for  a  portion  of  the  data 

presented  in  Chapters  5  and  6. 

Table  D.l :  Set-ups  for  measuring  temperature  profile.  Switch  thermocouple  readers  to  2nd 

for  2nd  half  of  experiment 


Location 


Zone  1  Exit 


Zone  2  Exit 


Zone  3  Exit 


Zone  4  Exit 


Zone  5  Exit 


Wall 


X 


X 


X 


X 


X 


.75 


A 


B 


C 


I) 


A 


B 


C 


D 


.25 

Center 

D 

C 

E 

D 

A 

E 

B 

A 

C 

B 

j  st  2nd 

Cool  down  and  Purge  with  C02  Between  Letter  set-ups 

Each  lettter  set  requires  4  runs  no  acoustics  and  with  acoustics  at  1 .0  LHSV-M  and  2.0 
LHSV-M 

Table  D.2:  Run  order  for  temperature  profile  (not  randomized). 


1 .0  LHSV-M 

2.0  LHSV-M 

w/acoustics 

w/o  acoustics 

w/acoustics 

w/o  acoustics 

Setup 

1st 

2nd 

,s, 

'jnd 

JS, 

^nd 

P 

^nd 

A 

1 

2 

J 

4 

B 

5 

6 

7 

8 

C 

() 

10 

11 

12 

I) 

13 

14 

15 

16 

E 

17 

18 

19 

20 

135 


136 


Table  D.3:  randomized  run  orders  for  catalyst  bed  length  and  parameter  settings  including 
replications. 


flow 


5  ml/min  setpoint 


10  ml/min  setpoint 


run 


Catalyst 

Bed 
Length 


full 


5/6 


2/3 


1/2 


1/3 


1/6 


w/  Acoustics 


10 


11 


11 


10 


2 
5 
8 
12 


w/o  Acoustics 

8 

12 

5 

4 

7 

10 

10 

1 

8 

9 

8 

4 

12 

4 

9 

8 

3 

6 

w/  Acoustics 

10 

6 

3 

12 

9 

3 

3 

12 

9 

6 

7 

2 

3 

6 

10 

9 

1 

5 

w/o  Acoustics 

7 

2 

9 

5 

6 

1 

11 

7 

4 

12 

11 

3 

2 

11 

7 

2 

7 

11 

Table  D.4  Run  Orders  for  Checking  Curvature  on  Flow  Rates. 


w/  Acoustics 

w/o  Acoustics 

flow  - 
(ml/min 
setpoint) 

6 

7 

8 

9 

6 

7 

8 

9 

run  -» 

4 

5 

8 

1 

3 

6 

7 

2 

Table  D.4  Run  Orders  for  Checking  Curvature  on  Sounc 

Pressure 

Pressure  setpoint 

0 

1000 

2000 

3000 

4000 

5000 

run  - 

4 

6 

2 

3 

5 

1 

APPENDIX  E 
TRANSFER  FUNCTIONS 


The  following  pages  are  the  transfer  functions  collected  from  the  function 
generator  to  the  pressure  transducer  signal  output  with  the  catalyst  bed  housing  only 
containing  33  °C  air  for  a  range  of  frequencies  of  40-2,100  Hz.  Similar  Transfer 
Functions  were  acquired  for  the  range  of  2,200-15,000  Hz. 


Function 
Magtfude 


Figure  E.1:  Transfer  Function  Magnitude  -  Pressure  Transducer  at  Port  1 


37 


138 


Figure  E.2:  Transfer  Function  Phase  -  Pressure  Transducer  at  Port  1 


Figure  E.3:  Coherence  -  Pressure  Transducer  at  Port  1 
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Function 
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Figure  E.4:  Transfer  Function  Magnitude  -  Pressure  Transducer  at  Port  2 


Tiarofef 
Function 
Phase  (deg)    .20.0-j 


2200.0  2400.C 


Figure  E.5:  Transfer  Function  Phase  -  Pressure  Transducer  at  Port  2 
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Coh»eiice 


00-, 
40.0         200.0 


400.0  600.0  800.0  1000.0  1200.0  14000  16000  1800.0  2000.0  22000  2400.C 

Hz 


Figure  E.6:  Coherence  -  Pressure  Transducer  at  Port  2 
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Figure  E.7:  Transfer  Function  Magnitude  -  Pressure  Transducer  at  Port  3 
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Ttarafa 
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Ph*re(dea]    .jgg 
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Hz 

Figure  E.8:  Transfer  Function  Phase  -  Pressure  Transducer  at  Port  3 


Figure  E.9:  Coherence  -  Pressure  Transducer  at  Port  3 


Tiarefei 
Function 
Magilude 


Figure  E.10:  Transfer  Function  Magnitude  -  Pressure  Transducer  at  Port  4 
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Figure  E.  1 1 :  Transfer  Function  Phase  -  Pressure  Transducer  at  Port  4 
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Figure  E.12:  Coherence  -  Pressure  Transducer  at  Port  4 
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Figure  E.13:  Transfer  Function  Magnitude  -  Pressure  Transducer  at  Port  5 
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Figure  E.14:  Transfer  Function  Phase  -  Pressure  Transducer  at  Port  5 
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Figure  E.15:  Coherence  -  Pressure  Transducer  at  Port  5 
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Figure  E.16:  Transfer  Function  Magnitude  -  Pressure  Transducer  at  Port  6 
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Figure  K.17:  Transfer  Function  Phase  -  Pressure  Transducer  at  Port  6 


146 


Coheience 


1.0- 


0.3- 


as- 


0.7- 


06- 


05- 


04- 


03- 


02- 


01- 


0.0- 


400         2O0.0  400.0  G0O0  8000 


1200.0  1400.0  1600  0         1800.0         2000.0  22000  2400  C 


Figure  E.18:  Coherence  -  Pressure  Transducer  at  Port  6 


APPENDIX  F 
TEMPERATURE  PROFILES 


The  following  pages  show  the  actual  averages  of  the  raw  data  points  collected  for 
temperature  and  the  regressions  used  to  construct  Figures  6.1  -  6.4.  Each  data  point 
shown  was  an  average  of  200  (minimum)  individual  data  points  collected  at  1  sample  per 
second.  This  appendix  presents  the  1 .0  LHSV-M  data  followed  by  the  2.0  LHSV-M  data. 
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